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ABSTRACT 


In this thesis, a scheme for the identification and classification of orthogonal 
frequency division multiplexing based signals is proposed. Specifically, the 
cyclostationary signature of IEEE 802.11 and IEEE 802.16 standard compliant 
waveforms is investigated. A model is introduced that identifies the waveform and, in 
the case of IEEE 802.11, confirms identification decision via cyclostationary feature 
extraction. If the waveform is identified as being IEEE 802.16 compliant, the scheme 
will classify the cyclic prefix size of the waveform. After cyclic prefix classification, the 
802.16 waveform is subjected to cyclostationary feature extraction for identification 
confirmation. The cyclostationary signature of each waveform is generated via a 
computationally efficient algorithm called the fast Fourier transform accumulation 
method, which produces an estimate of the waveform’s spectral correlation density 


function. Simulation results based on MATLAB implementation are presented. 
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EXECUTIVE SUMMARY 


Developing countries around the world are providing broadband wireless access 
to significant portions of their populations through Wireless Local Area Network 
(WLAN) and Wireless Metropolitan Area Network (WMAN) technologies. Many of 
these developing countries currently implementing broadband wireless networks are 
located in regions that possess potential or known threats to the security of the United 
States of America. This has lead to an increased interest in the ability to identify and 


classify these types of signals. 


There are three well-known methods of signal identification and classification 
currently being investigated for feasibility in this endeavor: energy detection, matched 
filter detection and cyclostationary feature detection. This thesis employs the 
cyclostationary feature method, due to the cyclostationary characteristics present in 
Orthogonal Frequency Division Multiplexing (OFDM) modulated signals. OFDM is the 
modulation scheme utilized in the IEEE 802.11-2007 WLAN and 802.16-2004 WMAN 


standards because of its resistance to multipath fading. 


The objective of this thesis is to develop a method of identifying and classifying 
OFDM based wireless data and communication network waveforms. In order to 
accomplish this, a scheme is developed to integrate the processes of initial waveform 
identification, cyclic prefix classification and cyclostationary feature extraction. 
Waveform identification is accomplished through a preamble correlation process. The 
results of the preamble correlation process determine if a received waveform is IEEE 


802.11 or 802.16 standard compliant. 


Cyclostationary feature extraction is employed to classify the cyclic prefix of 
IEEE 802.16 waveforms, as well as confirm the identification results of preamble 
correlation for both waveforms. The method of analyzing the cyclic spectral properties 
of OFDM signals is implemented through the fast Fourier transform Accumulation 
Method (FAM). The FAM algorithm produces an estimate of the analyzed signal’s 
cyclic spectral features. Additionally, the FAM provides a method of examining both the 


Xlll 


power spectral density and the cyclic spectral density commonly referred to as the 
spectral correlation density function. All results were generated through MATLAB 


simulation. 


There are four major contributions made by this thesis. First, the work 
established an approach to identify a received waveform with no a priori knowledge of its 
origin. After much research and experimentation, it was determined that a preamble 
cross-correlation operation provides satisfactory results while minimizing computational 


complexity of the simulations. 


Next, a method of classifying the cyclic prefix length of IEEE 802.16 waveforms 
is devised. Since the cyclic prefix option is established when the network is first setup, 
there is no way for a passive listener to identify this value by decoding captured control 
transmissions. By employing cyclostationary feature extraction, we are able to determine 


the cyclic prefix length in a rapid fashion. 


After identifying three methods of implementing cyclostationary feature 
extraction from [1], [3] and [4], we needed to determine which method or methods would 
be most compatible with the stated objectives of this work. Although embedded 
cyclostationary signature extraction is a promising technique, pilot subcarrier 


cyclostationary feature extraction proved to be the most applicable method. 
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I. INTRODUCTION 


A. BACKGROUND 


The rapidly advancing technologies of wireless communication networks are 
providing enormous opportunities. A large number of users in emerging markets 
throughout the world are gaining access to these networks. Developed countries possess 
extensive wired communications networks, but developing countries do not have the 
financial resources to implement the infrastructure necessary to provide access for 
significant portions of their populations. With the implementation of relatively 
affordable technologies based on the Institute of Electrical and Electronics Engineers 
(IEEE) Wireless Local Area Network (WLAN) and Wireless Metropolitan Area Network 
(WMAN) standards, these developing countries now have the ability to provide that 
access. Many of the developing countries currently implementing these broadband 
wireless networks are located in regions that possess potential or known threats to the 
security of the United States of America. This has lead to an increasing amount of 


interest in the ability to identify and classify these types of signals. 


There are three well-known methods of signal identification and classification 
currently being investigated for feasibility in this endeavor: energy detection, matched 
filter detection, and cyclostationary feature detection. This thesis employs the 
cyclostationary feature method, due to the cyclostationary characteristics present in 
Orthogonal Frequency Division Multiplexing (OFDM) modulated signals [1]. OFDM is 
utilized in the IEEE 802.11-2007 and 802.16-2004 standards. The communications 
industries refer to the implementation of these standards as Wi-Fi (IEEE 802.11) and 
Worldwide interoperability for Microwave Access (IEEE 802.16), respectively. Any 
further mention of IEEE 802.11 and IEEE 802.16 in this thesis will be in reference to the 
IEEE 802.11 and 802.16 standards, respectively. 


Traditional signal analysis methods typically employ stationary probabilistic 
models that examine a small time sample of a signal and assume that the statistical 


properties of the signals are time invariant. The term that identifies this stationary 


property with respect to first and second order statistical moments is Wide Sense 
Stationarity (WSS). The WSS property, which is satisfied if the signal’s mean is constant 
and autocorrelation is independent of time, ignores the statistically periodic or “cyclic” 
properties of manmade signals that can be developed. These cyclostationary properties 


can be very useful in distinguishing a variety of random signals from one another. 


B. OBJECTIVE AND APPROACH 


The objective of this thesis is to develop a method of identifying and classifying 
OFDM based wireless data and communication network waveforms. Waveform 
identification will be accomplished through a preamble correlation process. The results 
of the preamble correlation process will determine if a received waveform is IEEE 802.11 


or 802.16 standard compliant. 


Cyclostationary feature detection will be employed to classify the Cyclic Prefix 
(CP) of IEEE 802.16 waveforms and then confirm the identification results of preamble 
correlation. The method of analyzing the cyclic spectral properties of OFDM signals will 
be implemented through the Fast Fourier Transform (FFT) Accumulation Method 
(FAM). The FAM algorithm produces an estimate of the analyzed signal’s cyclic 
spectral features through MATLAB simulation. Additionally, the FAM provides a 
method of examining both the Power Spectral Density (PSD) and the cyclic spectral 
density commonly referred to as the Spectral Correlation Density (SCD) function. All 
results will be generated through MATLAB simulation. 


C. RELATED WORK 


Most of the background information that exists on cyclostationarity has been 
influenced by W. A. Gardner’s work. In particular, [2] is a comprehensive compilation 


of works on cyclostationarity. 


Although there is a significant amount of literature on cyclostationarity, there are 
only a few articles that extend this method to OFDM based signals. A significant amount 
of insight into properties of OFDM that might exhibit recognizable cyclostationary 


signatures are provided by [1], [3], and [4]... In particular, [1] focuses on a 
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cyclostationary signature produced by pilot subcarriers which are embedded within each 
OFDM symbol. Unlike the approach utilized in [1], which analyzed waveforms that 
applied a 2.5 dB gain to the pilot subcarriers, the waveforms investigated in this work did 
not amplify the pilot subcarriers instead an averaging process was implemented. We will 
investigate two additional cyclostationary signature approaches mentioned [3] and [4], 
although their application was deemed incompatible with the stated objectives of this 


thesis. 


D. ORGANIZATION OF CHAPTERS 


Chapter II provides a basic review of OFDM and cyclostationarity. The review of 
OFDM is by no means comprehensive and mainly focuses on the properties that are 
exploited for identification and classification purposes. The chapter will examine why 
OFDM modulation is becoming so popular and how it is utilized within the IEEE 802.11 


and 802.16 standards. Concluding the chapter will be a discussion on cyclostationarity. 


Chapter HI introduces a method by which IEEE 802.11 and 802.16 waveforms 
will be identified and classified. A model is proposed that will identify and classify IEEE 
802.11 and 802.16 standard compliant waveforms. Included in this discussion is a brief 
look at the FAM estimator that is employed to generate the SCD function of OFDM 
based waveforms. Chapter IV introduces the MATLAB implementation of the model 
described in Chapter HI. An overview of the program utilized to generate, identify and 
classify IEEE 802.11 and 802.16 compliant waveforms is performed. Included is a 
discussion of how the FAM is modified and improved to facilitate identification and 
classification of these waveforms. Results of MATLAB simulation of the proposed 


identification and classification scheme are presented. 


Chapter V will provide a brief summary of work performed and discuss 
significant results obtained and opportunities for future work related to identification and 
classification of broadband wireless networking waveforms utilizing cyclostationary 


techniques. 


The Appendix will consist of MATLAB code that was developed for simulation. 
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I. REVIEW OF OFDM AND CYCLOSTATIONARITY 


The most common metric utilized to measure the performance of communication 
networks is transmission bit rate. As the shear volume of data that is required to be 
transmitted over wireless channels consistently increases due to the development of more 
complex applications, the demand for higher transmission bit rates appears to be 
insatiable. Unfortunately, in a wireless channel, as the transmission bit rate increases, a 
phenomenon known as multipath fading increases as well, thus driving bit transmission 
errors toward unacceptable levels. The most effective method in combating multi-path 
fading at high bit transmission rates is multi-carrier modulation, more specifically 
OFDM. In the following sections, we will take a detailed look at the fundamentals of 
OFDM and how it is implemented by the IEEE 802.11 and 802.16 standards. Finally, the 


concept of cyclostationarity will be introduced. 


A. OFDM PHYSICAL LAYER FUNDAMENTALS 


The concept of multi-carrier modulation has existed since the early 1960s. It is a 
simple yet effective solution to the multipath fading imposed bit rate limitation of the 
channel. The effects of multi-path fading are negligible at relatively low bit rates. At 
higher bit rates, the period of a data symbol is reduced and becomes comparable to the 


duration of the delay spread in the channel leading to Inter-Symbol Interference (ISI). 


In order to overcome ISI, multi-carrier modulation breaks up the high bit rate 
serial stream that is to be transmitted into lower rate parallel bit streams. This is 
accomplished by dividing the serial bit stream of rate R into L parallel lower rate bit 
streams, of R/L, which are then modulated onto L subcarriers distributed within a 
specified frequency band. Once the L bit streams are modulated to their respective 
subcarriers, the modulated subcarriers are summed and transmitted. Figure 1 depicts the 


concept of multicarrier modulation. 


R bps 





Cost2nyeHL-1Ap 
Figure 1. Basic Multicarrier Modulation (After [5]). 


Effectively, the original bandwidth W of the high rate bit stream R is segmented 


into L subbands of bandwidth W/L. The subcarriers are separated by an optimal distance 


within the frequency band, referred to as orthogonality, to avoid inter carrier interference. 


This is accomplished through the effective use of digital signal processing techniques 


such as FFT and Inverse FFT (IFFT) in the baseband, prior to RF modulation. To 


illustrate the orthogonal subcarrier separation, Figure 2 shows the frequency response of 


adjacent subbands ina FFT. To visualize orthogonality, notice that when the response of 


any one subband is at its maximum, the collection of spurious responses from all the 


remaining subbands is zero. 


Frequency Responce 























0.1 O41 


a 
Normalized Frequency 


Figure 2. Orthogonal Subcarrier Spacing through IFFT Implementation. 
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The process of OFDM modulation is depicted in Figure 3. Initially, the high rate 
serial bit stream R is converted into L parallel lower rate bit streams. Next each bit 
stream is segmented into groups of bits called symbols. The size of the grouping is 
dependent on channel conditions; the better the channel state the larger the number of bits 
used to form the symbol and vice-versa. The purpose of the grouping is to prepare the 
bits for baseband modulation or symbol mapping. The number of symbol levels M is 


defined as 
Miao (2.1) 


where & is the number of bits per symbol. The symbol period is defined as 

T =Kki,, (2.2) 
where 7; is the bit period. The baseband modulation methods employed in most OFDM 
systems consist of Binary Phase Shift Keying (BPSK), Quadrature Phase Shift Keying 
(QPSK), Quadrature Amplitude Modulation 16 (QAM-16), and Quadrature Amplitude 
Modulation 64 (QAM-64) for k =1,2,4,6, respectively. The output of the baseband 


modulation is data symbol X,, a complex valued term with in-phase and quadrature- 


phase (I-Q) components. Next the L data symbols are input to an IFFT operation, given 


as follows [6], 


Nipper /2-1 ; 
x[n] = > ae (2.3) 


k=-N ppp 2 


where N_,,; 18 the size of the IFFT and Af is the subcarrier spacing in Hz. 


Typically, not all of the IFFT inputs are utilized as data subcarriers. A certain 
number of subcarriers are used as pilot subcarriers for channel estimation. Additionally, 
numerous lower and upper end subcarriers are set to zero to reduce adjacent channel 
interference. The number of pilot and null or guard subcarriers depend on the size of the 
IFFT and the standard governing the communication system. The pilot subcarriers 
transmit a pseudo-random sequence that is known by the receiver. This allows for 
determination of channel conditions and therefore which baseband modulation technique 
to employ. The pseudo-random pilot sequence will be an important feature for system 


identification as will be demonstrated later in this thesis. 
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The inputs to the IFFT that are not nulls or pilot subcarriers are utilized as data 
subcarriers. For instance, the system in Figure 3 illustrates the example of a 16 point 
IFFT with ten data subcarriers, two pilot subcarriers and three guard subcarriers. The 
guard subcarriers are located in the two most negative frequency indexed subcarrier 
positions and the most positive frequency indexed subcarrier location. The 0" subcarrier 
null is referred to as the “DC” null and serves as the center point or zero frequency 


subcarrier. 


Parallel 
To Serial, 
Addition Pulse 
Of Shaping, ya 
. : Cyclic Digital to aH) 
Input bit Symbol Prefix Analog 


stream. Mapping Conversion 





Pi Data Rate 
R bps Bh xq x(n) 


Figure 3. OFDM Block Description for 16 point IFFT. 


Figure 4 depicts the spectral layout of the guard, pilot and data subcarriers of the 
OFDM transmitter depicted in Figure 3. Additionally, Figure 4 depicts the subcarrier 
frequency offset index, which is a reference for where the subcarrier is positioned within 


the waveforms spectral bandwidth. 


Pilot Subcarriers 
Guard Subcarriers Guard Subcarrier 
Na re 
Subcarrier index:-8-7 -6 -5 -4-3-2-1 012345 67 ESE 
Data Subearriers DC Me Data Subecarriers 
Figure 4. OFDM Frequency Subcarrier Description of a 16 Point IFFT 


System. 


Next, a Cyclic Prefix (CP) or guard interval is inserted into the OFDM symbol 
that was formed from the output samples of the IFFT. The CP is a portion of the higher 
index IFFT output samples. These samples are copied and appended to the OFDM 
symbol as the leading portion. The purpose of the CP is to prevent ISI among OFDM 
symbols during transmission through a multipath fading channel; CP size is dependent on 


channel conditions and the governing standard. 


Referencing the IEEE 802.16 standard [7], CP options consist of 1/4, 1/8, 1/16, 
and 1/32 of the IFFT size. If channel fading conditions are severe, the longest CP size 
will be selected which would be 1/4 length of the 256 point IFFT or 64 samples long. 
Figure 5 illustrates how a 1/4 length CP on a 256 sample OFDM symbol is created. This 
is accomplished by taking the last 64 samples and appending them to the front of the 
OFDM symbol. Unfortunately, there is a price to be paid for the addition of the CP onto 
each OFDM symbol. On the receiver end of the transmission, the CP is discarded. 
Because the CP is redundant information that is discarded, it is pure overhead which adds 


no useful data to the transmission and reduces throughput. 


OFDM SYMBOL 


ler, _ Nerv -rle—_—_ Symbol Length= N,,- _,| 


4 


Figure 5. 1/4 Length CP Where the Last N,,,, /4 Samples are Copied and 
Appended to the Leading Edge of the Symbol. 






Now that the OFDM symbols have been formed and are ordered in a serial 
fashion, the samples are converted into a continuous-time signal by passing them through 
a digital to analog conversion and pulse shaping operation. The continuous time signal is 


then up converted to a radio frequency carrier and amplified for transmission. 


Since the fundamentals of the OFDM Physical Layer (PHY) that pertain to this 
thesis have been discussed, we will now take a very brief look at the IEEE 802.11 WLAN 


and 802.16 WMAN standards. The purpose of this discussion will be to focus on the 
portions of these two standards that will play a key role in the identification and 


classification of the respective OFDM signals. 


B. IEEE 802.11-2007 SPECIFICATIONS 


We have developed a basic understanding of OFDM in the previous section. The 


discussion will now focus on how it is employed within the IEEE 802.11 standard [8]. 


1. Physical Layer 


For the discussion of the IEEE 802.11 standard, we will begin with the physical 
layer (PHY). The IEEE WLAN standard divides its PHY into two functional sublayers: 
The Physical Layer Convergence Procedure (PLCP) and the Physical Medium Dependent 
(PMD). The PLCP allows the Medium Access Control (MAC) layer to operate with 
minimum reliance on the PMD. The PMD essentially performs the role of actually 
transmitting and receiving data. The PLCP organizes MAC Protocol Data Units (MPDU) 
in to a desired frame format. These frames are then handed off to the PMD for 


transmission. 


Figure 6 illustrates the structure of the PLCP frame. The PLCP preamble and the 
data OFDM symbols are the key points of interest of Figure 6, with respect to this thesis. 
These two items form the basis of an identification and classification scheme that will be 


introduced in Chapter III. 


| PLCP Header | 
[a $ $ $$ t=, 


RATE | Reserved] LENGTH Parity | Tail | SERVICE 7 Tail 
PSDU 





as Coded/OFDM | Coded/OFDM 
"~ | (BPSK, r= 1/2) | (RATE is indicated in SIGNAL) | 
ss; | 


"Simbel 
12 Symbols _|One OFDM Symbo Variable Number of OFDM Symbols 








Figure 6. PPDU Frame Format (From Ref [8]). 
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Each PLCP Protocol Data Unit (PPDU) frame begins with a preamble. Overall, 
the length of the preamble is equal to the duration of four OFDM symbols. Figure 7 
depicts the structure of the short and long preamble training sequences. The first half of 
the preamble is composed of ten short training sequences. The short preamble’s purpose 
is to provide for signal detection, Automatic Gain Control (AGC) convergence, diversity 
selection and timing acquisition to name a few of its contributions. The second half of 
the PLCP preamble consists of two “long training sequences” that aid in enhanced 
frequency acquisition for the receiver [8]. These two long sequences are preceded by a 
CP that is formed and appended in the same fashion as with a standard OFDM symbol. 
The first half of the preamble omits the CP. The short training sequences that compose 
the preamble are formed using a normalized sequence of QPSK I-Q data, given by 
S(k)=~13/6{0,0,1+ 7,0,0,0,-1—,0,0,0,1+ /,0,0,0,-1- j,0,0,0,-1-j 
,0,0,0,1+ 7,0,0,0,0,0,0,0,—-1—7,0,0,0,-1—7,0,0,0,1+ 7 (2.4) 
,0,0,0,1+ 7,0,0,0,1+ 7,0,0,0,1+ 7,0,0}, for —-26<k < 26. 
This sequence is processed by a 64-point IFFT to produce a short sequence of 16 time 


samples having duration of 0.8 ps: 


15(n)=ws(n) >> S (ke?) ; (2.5) 


where w, (t)is a window function. 


Generation of the long sequence is performed in a similar fashion but by using a 
sequence of BPSK data given by 
L(k)={,1,-1,-1,1,,-1,1,-L1,11,1,1,1,-L-111-11-11 
»1,1,1,0,1,-1,-1,1,1,-1,1,-1,1,-1,-1,-1,-1, (2.6) 
1,1,-1,-1,1,-1,1,-1,1,1,1,]}, for —26<k < 26. 











As with the short sequence, the long sequence is processed by a 64-point IFFT operation 


in the following manner 


a J2mkAf|n-T, 
r,(n)=w,(n) > L(k)el’ (7) (2.7) 
k=-26 
with 7, accounting for the 1.6-1s long CP located at the leading edge of the symbol. The 
duration of each long training sequence after the IFFT is 3.2 us. Figure 7 illustrates the 
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format and duration of the preamble structure for the 20-MHz channel spacing option. 
Additionally, Figure 7 shows the signal field and data OFDM symbols that follow. The 
signal field contains the data rate and length of the following transmission. Notice that 
after the nonstandard OFDM symbol format of the preamble, the subsequent data is 
encapsulated in standard OFDM symbols that are 4.0 us in duration. 


8+8=16ups 











10 x0.8=8 us x0.8+2x32=8.0ps 


fact Geld ee : ar ee 
} Fs, Fs, Fs, Fs, 18, 184%, Ts, hs, Tg Ge Pee 1 P 








4 4 ee 
Channel and Fine Frequency RATE SERVICE=DATA DATA 
Offset Estimation LENGTH 


+<—__—_—_> +> 
Signal Detect, Coarse Freq. 
AGC, Diversity Offset Estimation 
Selection Tuming Synchronize 


Figure 7. OFDM training Structure for 20 MHz Channel Spacing Option 
(After [8]). 


Now that there is a basic understanding of how the PLCP frames are constructed, 
a closer look at the construction of the OFDM symbols is warranted. The specific 
parameters of the three inter channel spacing options available in the 802.11 standard are 
listed in Table 1. Of note, the number of subcarriers is listed at 52. This number includes 
only the pilot and data subcarriers. Additionally, there are 11 guard subcarriers and one 


DC null, which results in a 64-point IFFT forming the basis of each OFDM symbol. 


An important feature of the OFDM symbol with respect to this research is the 
pilot subcarriers. They are positioned at frequency indexed subcarriers -21, -7, 7 and 21 
as shown in this sequence: 
P(k)={0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0, (2.8) 
0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,—1,0,0,0,0,0}, for —-26<k < 26, 
which excludes the six negative and five positive guard subcarrier indices but does 
include the DC null. The non-zero terms in this equation are the permanent values of the 
pilot subcarriers that are transmitted in each OFDM symbol created by an IEEE 802.11 


transmitter. 
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Information data rate 


Modulation 


6, 9, 12, 18, 24, 36, 48, 
and 54 Mb/s 
(6, 12, and 24 Mb/s are 
mandatory) 
(20 MHz channel 
spacing) 


BPSK OFDM 
QPSK OFDM 
16-QAM OFDM 
64-QAM OFDM 


3, 4.5, 6, 9, 12, 18, 24, 
and 27 Mb/s 
(3, 6, and 12 Mb/s are 
mandatory) 
(10 MHz channel 
spacing) 


BPSK OFDM 
QPSK OFDM 
16-QAM OFDM 
64-QAM OFDM 


(1.5, 3, and 6 Mb/s are 
mandatory) 
(5 MHz channel 
spacing) 


BPSK OFDM 
QPSK OFDM 
16-QAM OFDM 
64-QAM OFDM 


= 7 (64 states) K= 
convolutional code 


7 (64 states) 
convolutional code 


Error correcting code 4 states) 
conv mdi code 


ee 


Occupied bandwidth | 16.6 MHz 4.15 MHz 


Table 1: 





IEEE 802.11 Major OFDM PHY Parameters (From [8]). 


The pseudo-random property of the pilot subcarriers is produced by varying the 
polarity of the non zero values listed in Equation 2.8. This is accomplished by 


multiplying each of the permanent P(k) values by a pseudo random sequence. The 


following equation depicts the 127 element sequence that controls the polarity of the four 


pilot subcarriers 


n, (2) =€,1,1,1, -1,-1,-1,1, -1,-1,-1,-1, 1,1,-1,1, -1,-1,1,1, -1,1,1,-1, 1,1,1,1, 1,1,-1,1, 
1,1,-1,1, 1,-1,-1,1, 1,1,-1,1, -1,-1,-1,1, -1,1,-1,-1, 1,-1,-1,1, 1,1,1,1, -1,-1,1,1, 
-1,-1,1,-1, 1,-1,1,1, -1,-1,-1,1, 1,-1,-1,-1, -1,1,-1,-1, 1,-1,1,1, 1,1,-1,1, -1,1,-1,1, 
-1,-1,-1,-1, -1,1,-1,1, 1,-1,1,-1, 1,1,1,-1, -1,1,-1,-1, -1,1,1,1, -1,-1,-1,-1, -1,-1,-1}. 


(2.9) 


The signal symbol is the first OFDM symbol in the PLCP frame to include pilot 


subcarriers. In a Kronecker product-like operation, each pilot subcarrier in P(k) is 


multiplied by the first element of n,(/) and processed by the IFFT to produce the signal 


field symbol. The next symbol will include the second element of n,(/) being 


multiplied with the four pilot subcarrier values. This pattern will continue until 127 
OFDM symbols have been created. On the 128" OFDM symbol, the sequence n »(/) 
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will repeat as a cyclic extension for the remaining symbols of the frame. The cyclic reuse 


of n, () creates a pseudo random sequence for each pilot subcarrier which aids in 


channel estimation. 


2. IEEE 802.11 Medium Access 


Medium access in IEEE 802.11 is accomplished using a contention-based 
approach called Carrier Sense Multiple Access with Collision Avoidance (CSMA/CA). 
In CSMA, users or stations will listen or “sense” the medium to determine its state (busy 
or idle) prior to transmitting. The state of the medium and the rules governing the COMA 
algorithm will determine whether the station transmits or delays transmission. The 
foundation of this protocol is formed by partitioning time into slot times. All actions are 
initiated at the beginning of one of these slots. When a station has data to transmit, it will 
sense the medium for traffic. If the medium is idle, the station will delay an Inter Frame 
Slot (IFS) period of time, sense again and transmit if the medium is idle. If the medium 
is busy, the station desiring to transmit will wait until the medium is idle, delay one IFS, 
and then enter a random backoff phase before transmitting. The purpose of random 
backoff phase is to avoid simultaneous transmissions by waiting stations. If a 
transmission is sent while a station is in backoff mode, the delaying station will pause its 


backoff timer until the transmission is complete then restart the timer. 

In the above discussion, we presented some important features of the IEEE 802.11 
WLAN standard. We will now review the IEEE 802.16 standard. 
c. IEEE 802.16-2004 SPECIFICATIONS 


The description of the IEEE 802.16 specifications [7] will begin with the PHY 
and briefly cover the MAC functions that are of significance to this thesis. Only 
downlink specifications of the standard will be addressed as the uplink transmissions are 


not utilized for the identification and classification purposes in this thesis. 
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1; IEEE 802.16 Physical Layer 


IEEE 802.16 transmissions are structured in a time division multiple access 
(TDMA) frame based format. Figure 8 illustrates the standard Time Division Duplex 
(TDD) frame format for an 802.16 transmission. The base element of the 802.16 frame is 


the physical slot, having the duration 


(2.10) 


where f, is the sampling frequency. The number of physical slots per frame is 


dependent on the data symbol rate and the frame length. A 10-ms frame, for instance, is 
comprised of 6667 physical slots, which equates to 138 OFDM symbols for both the DL 
and the UL at a sample rate of 4 MHz. Figure 8 illustrates the frame structure of a TDD 
transmission format. TDD is the preferred method of two-way communications within 
the IEEE 802.16 standard and provides for a down link (DL) and an up link (UL). The 
DL transmission begins with the preamble followed by the Frame Control Header (FCH). 
From the FCH on, the rest of the OFDM symbols created are of the standard symbol 
format which will be explained in the following discussion. The FCH is similar to the 
signal field of an 802.11 transmission in that it describes the transmission rates and 
format to the subscriber stations. Next, the DL map and the UL map describe when each 
subscriber station will listen for their transmissions on the DL and when to transmit their 


data on the UL. 
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Figure 8. TDD Frame Structure (After [7]). 


Each frame begins with a preamble. The preamble is comprised of two OFDM 
symbols that are constructed from a sequence of QPSK I-Q data as follows [7] 





X,,(k) = {1-j, 1-j, -1-j, 1+), 1-5, 1-j, -1+4, 1-j, 1-j, 1-j, 14), -1-j, 14), 14), -1-j, 144, -14, -14, 1-j, 

-1+j, 1-j, 1-j, -14, 14), 14, 14, -14), 14, 14, 14, 14), -14, 149, 14), -1-j, 144, -14, -1-, 14), 

-1+j, 1-j, 1-j, -1-4, 143, 1-j, 14, -14), 14, 14, 14, 144, -1-j, 144, 14), -1-j, 144, -14), -14, 14, 
-14+j, 14), 14, 1-j, -14j, 144, 145, -1-j, 14), 144, 145, -145, 1-7), -145, -145, 13, -14, 1-j, 14, (2.11) 
1+, -1-j, -1-j, -1-j, -1+j, 1-j, -14, -1-j, 14, -1-j, -1-j, -1-j, 1-j, -1+), 1-j, 14, -14), 14, -1+4, 
-14), -1-+j, 14, 0, -1-j, 1+), -149, -14), -1-j, 14), 1+), 145, -1-j, 14), 1-j, 1-j, 1-j, -14j, -1+j, 
-14+j, -14j, 1-j, -1-j, -14, -14), 1, 14), 144, -14+), 14, 14, 1-j, -14j, 14, -14, -14, -13, 144, 
1+j, 14, 144, -14, -14), -14j, 144, -13, 14, 14, 145, -1-j, -1-, -1-j, 14), -14, -14, -144, 
-14j, 1-j, 1-j, 1-j, 1-j, -1+j, 14), 14), -1-j, 144, -14, -14, -1-j, 14j, 144, 14, -1-j, 14), 15, 
1-j, 1-j, -1+, -14), -144, -14j, 1-j, -14, -1-j, 1-j, -149, -15, -1-j, 15, -14, -14), -144, 14, -14), 
1+j, 14), 14), -14, -14, -14, -14, 145, 14, 1-3} 
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The first preamble symbol consists of four short sequences generated by 


decimating YX, (x) according to the following equation 


s(k)= 2X7,(k) kmod4=0 (2.12) 
0 kmod4 +0. 


After the decimation process is performed, the 50 frequency samples generated are 
concatenated with three other identical groups of 50 samples. This group of 200 samples 
is now subjected to the IFFT operation with the lower, upper and DC nulls inserted into 


the appropriate locations. Figure 9 illustrates this process by showing the format of the 
16 


samples prior to the IFFT and then the OFDM symbols they form after the IFFT 
operation. | The second preamble symbol is created in a similar fashion. Instead of 
creating 4 short sequences, two long sequences are produced by the following modulo 2 


decimation process: 


2X p(k Z= 
L(k)= V2X p(k) kmod2=0 (2.13) 
0 kmod2 #0. 
These two long sequences are 100 samples in length and are displayed in Figure 9. 
Nulled Subcarriers Nulled Subcarriers 


™ ea 









Figure 9. Downlink and Network Entry Preamble Structure (After [6]). 


Following the preamble symbols within an IEEE 802.16 frame are the Frame 
Control Header (FCH), the DL MAP and UL MAP. These control messages are 
embedded in OFDM symbols that are of the same format as the data carrying symbols. 
The construction of the data carrying symbols is a key feature of our waveform 


identification and classification method and will be the next topic of discussion. 


An IEEE 802.16 OFDM symbol is constructed in a fashion similar to that of 
802.11. The dimensions are the main difference. For instance, the lower and upper 
guard nulls consist of 28 and 27 subcarriers, respectively. There are eight pilot and 192 
data subcarriers. Including the DC null and there are 256 subcarriers in an IEEE 802.16 
OFDM symbol. Unlike IEEE 802.11 waveforms, the 802.16 OFDM symbol CP length is 
not a fixed size. In addition to the 1/4 length CP, which is standard in 802.11, there is the 
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option of a 1/8, 1/16 and 1/32 length CP with the 802.16 standard. These shortened CP 
options reduce excessive overhead in the event multi-path fading is not severe enough to 


justify the a 1/4 CP. 


The pilot subcarriers are modulated with a pseudo random binary sequence 
(PRBS) that is created by the eleven-bit shift register depicted in Figure 10. The output 
of this PRBS generator is show by the following binary irreducible primitive polynomial 
[7], 

w(k)=X" 4X? +1. 


The final values that are loaded into the pilot subcarriers prior to the IFFT are determined 
in the following fashion, 


Pog (k) = Psg(K) = P,(k) =P, (k)=(1-2w(k)) 


Pak) =Pas(k) =P (k) = Bs (Kk) =(1-20(4)) (2.14) 


where the subscript of P(k) represents the corresponding pilot subcarrier frequency 


offset index and w(k) is the compliment of w(k). The pilot subcarriers are loaded with 


the results of Equation 2.14, OFDM symbol by OFDM symbol. 





Figure 10. PRBS generator for IEEE 802.16 Pilot Subcarrier Modulation 
(After [7]). 
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The baseband modulation schemes used in IEEE 802.16 waveforms are the same 


four options discussed in Section B.1. 


Now, we will take a brief look at the MAC scheme employed by the IEEE 802.16 


standard. 


2. IEEE 802.16 Medium Access 


The MAC scheme utilized by the IEEE 802.16 is Time Division Multiple Access 
(TDMA), a contention-free approach. This is a very structured scheme where all users 
communicate over the same channel during an assigned window in time. A centralized 
controller (the base station) establishes a schedule identifying when each user will listen 
for their traffic on the DL and when they will transmit on the UL. With low user volume, 
this method can be quite wasteful. However, as the number of users increases, this 


method becomes much more efficient than contention based approaches. 


Figure 8 depicts the structured nature of a TDMA-based medium access scheme. 
Each frame begins with a DL transmission that is roughly half of the frame length. 
Following the preamble, the FCH informs the Subscriber Stations what baseband 
modulation and coding is being used on their transmission. The DL MAP tells the 
Subscriber Stations when it is their turn to listen for their data. The UL MAP informs the 


SS when it is their turn to transmit. 


At the beginning of each UL subframe, there is a contention window for initial 
ranging and new Subscriber Stations network entry. These bursts are short in duration 
and only consist of a few OFDM symbols. For this reason, we decided to focus on the 


more structured DL subframe for identifying and classifying 802.16 based waveforms. 


The discussion will now cover the fundamentals of cyclostationarity. 


D. CYCLOSTATIONARITY 


Traditionally, analysis of communication signals embedded in noise has assumed 
the statistical property of stationarity. Stationarity has several definitions but the most 
common is Wide Sense Stationarity (WSS). A random process (noisy communication 


signal) is said to be WSS if its first moment (mean) is constant and its second moment is 
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a function of the time difference. Even if these two conditions do not hold true for the 
entire time duration of the signal, small segments of the signal can have a constant mean 
and time invariant second moment [9]. Cyclostationarity, on the other hand, states that 
most manmade signals such as communication, radar and sonar signals, to name a few, 
possess first, second and up to nth-order statistical parameters that vary with time [2]. 
These variations can be periodic which implies cyclostationarity, or they can possess 
many different periods, which implies polycyclostationarity. The time varying statistical 
parameters can manifest themselves as additive sine wave components in the time 
domain or spectral lines in the frequency domain. Whether these periodicities exist in a 
cyclostationary form or a polycyclostationary form, they can provide an outstanding 


method for detecting and identifying different types of communication signals. 


Before the discussion goes into how these cyclostationary properties are 
exploited, a brief review of pertinent terms and properties will be performed. Statistical 
analysis of random processes is normally performed from two differing approaches: 
ensemble averages and time averages. The ensemble average approach involves analysis 
of a sample space of repeated events. Generally, this utilizes density and distribution 
functions that characterize the ensemble. The time average approach, often referred to as 
the time-series approach, deals with averages over the time axis of a random signal. The 


ensemble average is given by 
E{g[X()}4 | [XO] (xa, (2.15) 
where E {+} is the expectation operator and g(-) is a function. The time average 


operation is given by 


Egle} Stim Lf elx(e+r)y er, (2.16) 


where E {| the time averaging operator. These two statistical analysis approaches can 


be equated through the principal of ergodicity. If a random process is ergodic, the 


statistical moments generated via ensemble and time averages are equivalent. 
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In order to develop the cyclic moments that identify a signal’s cyclostationarity, a 
quadratic transformation is performed. This operation is not unlike the autocorrelation 


function described as follows 
R,(t,t-1) =E{x(t)x’ (¢-7)}. (2.17) 
To develop the fundamental parameter of second order cyclostationarity, the Cyclic 


Autocorrelation Function (CAF), the process in Equation 2.17 is performed in 


conjunction with the Fourier Transformation of x(¢) in the following fashion 
3 A A y - jano{ 1-5] 
ta (=z x(t)x (t-r)e : (2.18) 


where a is referred to as cyclic frequency and is displayed on a separate axis. This 
separate axis allows for the visualization of the amount of correlation that exists between 
frequency shifted versions of the signal. When a is equal to zero, the CAF becomes a 
traditional autocorrelation function. As a becomes non zero, and cyclic properties 


possessed by x(t) will be displayed along the cyclic frequency axis. 
To better illustrate this cyclic phenomenon, consider an amplitude modulated 


signal [2] 


x(t)= x4 (t) cos(27 ft + 0) 
i (2.19) 
eco i(2af,t+0) -j(27f,t+9) 
=5a(s)[e’ +e! | 


where a(t) is a random signal whose autocorrelation is nonzero and possesses the 


following property 
E{a(t)a'(t-r)e™\=0 forall a #0. (2.20) 


The next step will be to calculate CAF of x(#) as follows 
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(2x f,t) o * —j2aat 1s (2x f,t) ~ * —j2nat 
e! E{a(t)a (t-r)e? eae E{a(t)a (t-r)e” \ 
ef Fe iho Ea (t) a (t — T) era (2.21) 


e179 gi?*hol E {a (t) a (t = T) e 12ala+2 fo) } ; 


Applying the property of Equation 2.20 to Equation 2.21 yields 
1 


re (z) for a =+2f, 
Re (¢)= $R, (z)o0s (22 f,2) for a =0 (2:22) 
0 otherwise. 


Notice, whena@ =0, Equation 2.22 is equal the stationary autocorrelation of x(t). Now, 
we will extend this analysis to the frequency domain to observe the spectral properties of 
x(t) : 


According to Gardner [2], by taking the Fourier transform of a continuous time 


CAF, we produce a function called the Spectral Correlation Density (SCD) as follows 


Sef \= y Reger, (2.23) 
By applying Equation 2.22 to oe 2.23, we have 
weirs, (/) for a =+2f 
S?(f)= TSF +S)4+ 25 (F-f,) for a =0 (2.24) 
0 otherwise. 


Notice that when a =0, the SCD of x(t) reduces to the PSD for this signal. 


This chapter covered the fundamentals of OFDM modulation. It described how 
the IEEE 802.11 and 802.16 standards utilized OFDM modulation to produce 


communication waveforms and transmit them over a shared wireless channel. Finally, 
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the concept of cyclostationarity was introduced. The topics discussed in this chapter will 
be built upon in the following chapters, as we discuss the methods utilized to identify and 


classify OFDM based wireless network waveforms. 
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HiIl. OFDM BASED SIGNAL IDENTIFICATION AND 
CLASSIFICATION 


Currently, most modern wireless data communication networks are utilizing 
OFDM based modulation techniques of differing configurations. With this in mind, a 
reliable waveform identification and classification method is becoming necessary. It is 


especially imperative from an intelligence collection perspective. 


Within in the field of cognitive radios, there is much research being conducted on 
the topic of waveform identification and classification. Cognitive radios maximize usage 
of idle portions of the electro-magnetic spectrum. When the licensed owner of a reserved 
frequency band is not actively using this resource, cognitive radios will sense the 
unutilized portion of spectrum and make use of it [4]. Once the licensed user begins to 
transmit, the cognitive radio must recognize this traffic and vacate that frequency band. 
Of the three waveform identification methods mentioned earlier, the trend appears to be 
toward cyclostationary feature detection. Although matched filter, and energy detection 
methods hold promise in waveform identification, the matched filter requires a very 
complex hardware footprint to cover all possible waveforms. Energy detection is 
susceptible to noise and interfering signals [1]. Cyclostationary feature identification 


provides for precise and reliable waveform identification and classification. 


In this chapter, we will introduce a model to identify and classify IEEE 802.11 
and 802.16 waveforms. A method of waveform identification by preamble correlation is 
discussed. Next, three methods of cyclostationary feature identification are presented and 
examined in detail. Finally, a computationally efficient algorithm is introduced that 


produces a SCD estimate of the waveforms analyzed in this work. 


A. WAVEFORM IDENTIFICATION AND CLASSIFICATION 


Waveform identification and classification can conjure different meanings for 
different applications. In the context of the thesis, identification shall refer to the 
determination of which wireless communication standard the waveform is compliant 


with. Classification shall refer to the determination of which variant of cyclic prefix (CP) 
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is utilized when constructing the OFDM symbols of the signal being analyzed. In order 
to achieve the stated objective of this thesis while complying with the above-mentioned 
definitions, a waveform identification and classification model is proposed as shown in 
Figure 11. The identification aspect here is limited to the IEEE 802.11 and 802.16-2004 


standards. 


Waveform 





Passive Signal Preamble 
Capture Correlation 





Figure 11. Proposed Waveform Identification and Classification Model. 


Once a signal of interest has been captured, it will be down converted to baseband 
for analysis. Since the signals analyzed in this work are compliant with published IEEE 
standards, the beginning of a transmission burst or frame will start with a preamble. The 
purpose of the preamble is to allow the receiver to detect the presence of a compliant 
signal and for synchronization in both time and frequency. In order to determine if the 
captured signal is an IEEE 802.11 or an 802.16 waveform, a cross-correlation operation 
is performed. If the received waveform is identified as 802.11, the identification decision 
will be confirmed by cyclostationary analysis leading to the classification of the signal. 
Since the 802.11 standard provides for only one CP length, CP classification is not 
warranted. When an 802.16 signal is identified, the first step in the signal classification 
process is to determine the CP of the waveform because of the potential of multiple CP 
sizes. In order to generate a cyclostationary signature of the waveform, the CP size must 
be determined and those samples removed prior to FAM estimator. After the CP size has 
been determined, the cyclostationary signature is extracted and the waveform is 


classified. 
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The following discussion will explain the process by which a received waveform 


is identified through a preamble cross-correlation operation. 


B. WAVEFORM IDENTIFICATION BY PREAMBLE CROSS- 
CORRELATION 


This section focuses on the method employed to identify IEEE 802.11 and 802.16 
standard compliant waveforms through a cross-correlation process. In this discussion, 
two signals or processes will be considered. One will be a reference signal x(t) and the 
other will be a received signal y(t) with additive white Gaussian noise (AWGN). In order 
to determine if these two signals possess some degree of statistical similarity, the time 


averaged cross-correlation function is calculated as follows [9], 


Rev (r)= EX" (t)¥(—2)} = him = x(0)y(e—e)de. G.1) 


This operation is referred to as the time averaged cross-correlation function. Assuming 
ergodicity, the corresponding ensemble cross-correlation function can be obtained 


through the expectation operation, 
Ry (t)=£{X"(t)¥(t-7)}. (3.2) 
For waveform identification purposes, let x(t) represent an IEEE 802.11 standard 
compliant preamble sequence consisting of 320 samples, s,(t), ic. x(¢)=5,(¢). Let 
y(t)be equal to the same IEEE 802.11 preamble sequence with AWGN, ic., 
y(t)=s,(t)+n(t) where n(t) represents noise and s,(t) and n(t) are uncorrelated. 

The cross-correlation of these two processes given by 

Re (t)= Elsi Os +a)}, G3) 


which reduces to 


Ry (t)=E{s; (t)s,(t—-7)}. (3.4) 
The cross-correlation function for zt = 0 is given by 
Re (0)=E ls (0 } (3.5) 
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In the case where x(t) is a 320 sample IEEE 802.16 preamble sequence, i.e., 
x(t)=5,(t)and s,(t)and n(t) are uncorrelated, Equation 3.2 reduces to 
Ry (0)=E{s; (t)s,(o)}. (3.6) 


In practice, the discrete time cross-correlation function based on Equation 3.1 is 


computed as follows 


N-1 


R, [== > x [Ap[A+ 7]. (3.7) 
N yn 
Following on the lines of the assumptions made in arriving at Equations 3.5 and 3.6, we 
have 
Toe ; 
Rr [0]=+ ¥ fo] G.) 
k=-N 
when x(t)=s,(t) and 
1 N-1 
Ry [0J=— & s3[n]s,[n] (3.9) 
N =n 


when x(t)=s,(t). 


Figure 12 shows a MATLAB generated result of the cross correlation between an 
IEEE 802.11 standard compliant preamble sequence versus a received signal from an 
AWGN channel and 802.16 standard compliant preamble sequence versus the same 
received signal from an AWGN channel. Each sample preamble is comprised of 320 
samples, which equates to the entire 802.11 preamble sequence (short and long training 
sequences). Only the first half of the 802.16 preamble (first four 64 sample sequences 
with 1/4 length CP) is included in a sample window of this size. In addition, the signal- 
to-noise ratio (SNR) for the received signal was selected to be 0 dB. Within the context 


of this work, SNR is defined as 


SNR =10log,, ele Power 
Noise Power 
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At /=0, Figure 12(a) demonstrates a very high degree of correlation. This indicates a 
match with the first 320 samples of the received signal. The result displayed in Figure 


12(b) indicates little to no correlation. 
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Figure 12. Cross-Correlation of an IEEE 802.11 Preamble Versus (a) 802.11 


Sample Preamble, (b) 802.16 Sample Preamble. 


Figure 13 depicts the results when an IEEE 802.16 waveform is received from an 
AWGN channel and subjected to the same cross-correlation process. As is the case with 
the results displayed in Figure 12, there is an identifiable match and a definitive 
mismatch. The 802.11 cross-correlation appears to display little to no commonality while 
the 802.16 preamble cross-correlation indicates a high degree of likeness. For a much 
more thorough discussion on the process of signal detection and classification via IEEE 


802.11 and 802.16 preamble cross-correlation refer to [12]. 
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Figure 13. Cross-Correlation of 802.11 Preamble Versus a) 802.11 Sample 
Preamble, b) 802.16 Sample Preamble. 





Now that the preamble cross-correlation operation has identified the type of 
waveform that was received, the next step will be to confirm the preamble identification. 
This will be accomplished via cyclostationary feature extraction. In the event an IEEE 
802.11 waveform was identified, the cyclostationary feature extraction will confirm the 
results of the preamble cross-correlation. This process will entail passing the received 
signal into the FAM SCD estimator and confirming its cyclostationary signature and will 
be discussed in the next section. If an IEEE 802.16 waveform is identified, 
cyclostationary feature extraction will be used to determine which CP length was used to 
construct the OFDM symbols. Once the CP length is determined, this signal will be 
processed by the FAM algorithm to confirm the waveform identification results of the 


preamble cross-correlation and the CP classification results. 


30 


C. CYCLOSTATIONARY FEATURE EXTRACTION 


As mentioned in Chapter I of this thesis, there are three widely known methods of 
cyclostationary feature extraction in reference to identification and classification of 
OFDM based signals. The first method exploits the repetitive transmission of the frame 
preamble [3], the second utilizes or takes advantage of the cyclostationary signature 
generated by the pilot subcarriers imbedded in each OFDM symbol [1], and the third 
embeds a predetermined cyclostationary signature into each OFDM symbol, which can 
be recognized by cyclostationary analysis [4]. In the following discussion, we take a look 


at each on of these options under consideration. 


if Frame Preamble Cyclostationary Signature 


As mentioned in the previous chapter in discussions on the IEEE 802.11 and 
802.16 frame construction, all burst transmissions of frame data are preceded by a 
preamble. The distinct pattern possessed by preambles from each standard allows a 
receiver to detect the presence of a standard compliant signal, rapidly synchronize to the 
signal and ultimately recover the data imbedded in the OFDM symbols. At first glance 
this would appear to be a solid method of waveform identification by cyclostationary 
analysis [3]. Unfortunately, each standard possesses structural elements that make this 


approach unsatisfactory. 


a. IEEE 802.11 Preamble Considerations 


In the IEEE 802.11 standard, frame transmissions do not occur in a set, 
periodic fashion. In fact, considering that frame lengths vary in size, especially when 
comparing the downlink and uplink transmissions, and the mechanics of CSMA/CA 
access scheme, preamble transmissions are essentially a random phenomenon within a 
WLAN network. This leads to the conclusion that there was no value added by this 


identification method over a conventional cross-correlation of preamble time samples. 


b. IEEE 802.16 Preamble Considerations 


The limitations with respect to the IEEE 802.16 standard have to do with 


the variable length CP appended to the beginning of each OFDM preamble symbol. The 
31 


CP tends to disrupt the cyclostationary signature calculation. Since the CP can be one of 
four possible values depending on channel conditions, as with the 802.11 standard, there 
is little gain to justify the increased processing required with generating a SCD of the 


preamble when the cross-correlation of time samples is sufficient. 


Ze Pilot Subcarrier Cyclostationary Signature 


The pilot subcarrier cyclostationary signature method focuses on the spectral lines 
formed in the cyclic frequencies by the pilot subcarriers. Of the three cyclostationary 
feature extraction methods investigated in this thesis, pilot subcarrier cyclostationary 
signature displayed the most promise. Pilot subcarriers are positioned throughout the 
bandwidth of each standard’s OFDM symbols in order to estimate channel conditions [6] 
[7]. The positioning of these pilot subcarriers is performed in a fashion to ensure the 
entire spectrum of each OFDM symbol is adequately assessed while minimizing 
overhead. If too many pilot subcarriers are dedicated to channel estimation, there are 
fewer data subcarriers dedicated to throughput. With these pilot subcarriers permanently 
positioned at the same subcarrier frequency offset index and consisting of a pseudo- 


random BPSK sequence, a cyclostationary pattern is formed. 


Ideally, the polarity of the pilot subcarrier BPSK sequence would be fixed. This 
would generate a robust cyclostationary signature at the pilot subcarrier positions of the 
principal cyclic frequency axis. Unfortunately, this would create a spectral signature that 
would manifest itself as spectral lines along the principal frequency axis of the signals’ 
PSD, which is undesirable as it would lead to a high peak-to-average ratio. A high peak- 
to-average ratio would become an issue when the signal is amplified prior to transmission 
[5]. For this reason, a pseudo-random BPSK sequence is modulated on to the pilot 
subcarriers. Unfortunately, the pseudo-random pilot subcarrier sequence requires the 
SCD of multiple OFDM symbols to be averaged in order to produce a strong cyclic 


signature for feature extraction. 


Figures 14(a) and 14(b) display the SCD and PSD of an IEEE 802.16 waveform, 
respectively. This result was generated by the FAM SCD estimator, which will be 


introduced in the next section. The SCD demonstrates a clearly identifiable pilot 
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subcarrier cyclostationary pattern. This is the signature that is exploited when classifying 
IEEE 802.11 and 802.16 waveforms, although it is not nearly as robust when the pilot 
subcarrier sequence is pseudo-random. Unfortunately, Figure 14(b) also exhibits spectral 
line development at the pilot subcarrier locations of the PSD. This will aggravate the 
peak-to-average ratio when the signal is passed through a non-linear amplifier prior to 


transmission. 
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Figure 14. Spectral Line Development: (a) 802.16 Waveform SCD with 
Constant Pilot Subcarrier Sequence. (b) 802.16 Waveform PSD with 
Constant Pilot Subcarrier Sequence. 


To reduce the peak-to-average ratio, both waveforms utilize a pseudo random 
BPSK sequence for transmission on the pilot subcarriers. The cyclostationary signature 


obtained is useful in identifying and classifying OFDM signals. 
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os Embedded Cyclostationary Signature 


The embedded cyclostationary signature method of OFDM based waveform 
identification is the topic of much research within the field of cognitive radios. The 
reason for all the interest in this method is its ability to transmit a unique identification 
signature. This signature could be employed by licensed users of a range of spectrum to 
alert non licensed users that the registered user is transmitting. By embedding a unique 
signature within the OFDM symbols that are carrying the data of the licensed user as 
underlying periodicities, normal communications can occur with little overhead in terms 


of bandwidth and power. 


This method is similar to the pilot subcarrier sequence method in that the receiver 
has a priori knowledge of a registered user’s signature, which is embedded within its 
waveform. In order to accomplish this, an operation referred to as subcarrier set mapping 
is performed [4]. Subcarrier set mapping is accomplished by mapping or routing the 
same I-Q baseband data inputs to “mirror” locations of the IFFT operation. The term 
“mirror” refers to the same subcarrier frequency offset indexes with respect to the DC 
null, i.e., data subcarriers 44 and—44. Subcarrier set mapping produces a cyclic spectral 
signature that is very distinct for different subcarrier set mapping schemes. Figure 15 is 
an example of subcarrier set mapping. This is one of many possible configurations of 


subcarrier set mapping. 


A major inefficiency of this scheme is illustrated in Fig. 15. In order to make the 
signature unique amongst a large population of licensed users, the number of mapped 
subcarriers must increase. The net effect of increasing the number of mapped subcarriers 
is the reduction of system throughput. This becomes highly undesirable as the number of 


mapped subcarriers begins to crowd out the number of data subcarriers. 
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Figure 15. An Example of Subcarrier Set Mapping to Establish an Embedded 
Cyclostationary Signature. 


Within the field of cognitive radios and spectral reuse schemes currently being 
explored, embedded cyclostationary signature provides a powerful method to positively 
identify all licensed users. Unfortunately, this method requires a cooperative population 
of users that will transmit their predetermined signature. This will most certainly not 
always be the case. In addition, the cyclostationary signature produced will allow for the 
identification of a specific transmitter (with a priori knowledge of the signature) but not 
the type of standard compliant waveform produced by the transmitter. With out this 


knowledge, recovery of the user data would be impossible. 


After weighing the advantages and disadvantages of each cyclostationary feature 
identification method, we determined that the pilot subcarrier cyclostationary signature 
method would be implemented in this work. The frame preamble cyclostationary 
signature method is excluded because the cross-correlation process is just as effective 
with significantly less processing requirements. The embedded cyclostationary signature 
method, while possessing potential in an identify friend or foe role, is not applicable to 
identifying transmissions from uncooperative parties; hence, it is not employed in this 


work. 


Next, the discussion will focus on the method utilized to produce the 


cyclostationary properties of an OFDM based signal. 
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D. FFT ACCUMULATION METHOD (FAM) 


As modern communication system waveforms have increased in complexity, a 
need arose for computationally efficient methods of producing the SCD of these signals, 
if their cyclic properties were to be exploited. An algorithm was developed in reference 
[10] that addressed this need. Two classes of cyclic spectral analysis algorithms were 
identified in the literature: frequency smoothing algorithms and time smoothing 
algorithms. While both classes of algorithms are effective at estimating a signal SCD, the 
time smoothing approach is considered to be more computationally efficient. Within the 
time smoothing class of algorithms, reference [10] developed two computationally 
efficient algorithms: the FAM and Strip Spectral Correlation Algorithm (SSCA) while 
reference [11] developed MATLAB code that implemented the FAM and SSCA 
algorithms for the detection and identification of cyclostationary signals. After reviewing 
the results in [11] and testing both algorithms by MATLAB simulation, it was determined 
that the FAM implementation of the time smoothing approach provided better results for 
the purpose of this work. This determination resulted from the better resolution provided 
by the FAM estimator. When considering the number of subcarriers involved in the 
construction of an IEEE 802.16 OFDM symbol, detailed resolution becomes an 


imperative feature. 


The basics of the time smoothing algorithms for computing the cyclic cross 


periodogram are described below. The cyclic cross periodogram is given by 


SZ (nf), =2E(X, (nf +a/2)¥ (nf —a/2)}, (3.10) 


where T is the length in seconds of a data tapering window which slides over a signal 


segment of duration At and X,(n,f+a/2) and Y,(n,f—a/2) are the spectral 
components of the sampled signals x(n) and y(n), respectively. These spectral 


components are referred to as complex demodulates and are generated by passing the 
signals through a narrow-band, band-pass filter. The complex demodulates are produced 


as follows 


N12 


X,(n,f)= > a(k)x(n-kye PU (3.11) 


k=-N12 
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and 


N'Y2 


Y) al(k)xe(n-k)e PO | (3.12) 


k=-N12 


Y,(n, f) 


where a(k) serves as a data tapering window of length T=N'T,, T, is the sample 


duration and N' is size of a sliding point FFT which is described in the next paragraph. 
Once computed, the complex demodulates are cross-correlated over a time span of 


At seconds. 

In an effort to reduce the computational complexity involved with the time 
averaging process, the FFT can be used to produce the SCD estimate of x(n) and y(n). 
Development of the complex demodulates remains relatively unchanged with the 
exception of a sliding N'-point FFT being incorporated to channelize the input samples. 


Channelization consists of blocks of Z samples being selected from the input data and 
transformed by the FFT in a decimation process intended to increase estimation 
efficiency. The value of N' is determined via the selected sampling frequency /, and 


Ss 


the desired frequency resolution Af in the following fashion 


f 

log,| <— 
yoxg! (a) (3.13) 
Next, the complex demodulates are downshifted in frequency to baseband. Now that the 
complex demodulates have been channelized and reside at baseband, the Y, (nL, f,) 
demodulate is conjugated and multiplied with X,(nZ, f,) in preparation for the P -point 


FFT. The value of the P is determined as follows 





f 

logy - 
ep) (real) (3.14) 
where Aq is the desired cyclic frequency resolution. Now substituting these parameters 
into Equations 3.11 and 3.12, the FAM estimator using FFT for implementing the time 


averaging process is given by [10] 


—j2arq 


Se" (nb. f,), = Xr (LAG (bf )g.(a-re ? Ct) 
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where g.(n-r) is the Hamming window operation, k and/=1,...,.N' and q is the 
channel index in the range 
PL re 
- <q< —1., 3.13 
2N' 2N' oe 


Figure 16 depicts a block diagram of a generic implementation of the FAM time 
smoothing algorithm. 
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Figure 16. Block Diagram FFT Accumulation Method. (From [11]). 


Since the process of generating an SCD estimate has been discussed, the 
discussion will now focus on how the results of the FAM estimator are presented. 
Normally, when presenting PSDs of a random signal, the standard method is to display 
the PSD as a function of frequency on a two dimensional graph. When displaying the 
SCD of a signal, an additional dimension must be included to represent the cyclic 
frequency axis. As described in [2], a region of support is formed within a bifrequency 
plane to represent a signal’s SCD function. Figure 17 is an illustration of the bi- 
frequency plane and the region of support for a lowpass signal. Note that the cyclic 
frequency bandwidth of the region of support is twice the conventional frequency 
bandwidth of a low pass signal. The term a represents the cyclic frequency axis and f 
represents the traditional frequency axis. In order to display the SCD produced by the 
FAM estimator, a region of support is populated within a N'+1x2N-+1 array, where 


N=PL is size of the data vector processed by the FAM. The dashed boundary 
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represents the dimensions of the array while the shaded area represents the cyclic region 
of support. Not all the cells of this array contain useful data, only the ones that fall 


within region of support. 





Figure 17. Region of support within the bifrequency plane (From [2]). 


The preceding discussion was intended to provide a basic understanding of how 
the SCD estimate is generated and displayed by MATLAB simulation. A more detailed 
discussion of the FAM algorithm will be conducted in Chapter IV. 


This chapter introduced a method for identifying and classifying OFDM based 
waveforms. A preamble cross-correlation operation was introduced to identify whether a 
waveform is compliant with the IEEE 802.11 or 802.16 standard. Next, three methods of 
cyclostationary feature extraction were introduced. Pilot subcarrier cyclostationary 


signature identification was determined to be the most applicable method of the three in 
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achieving the stated goal of this work. Finally, a computationally efficient method of 


producing the SCD of the OFDM waveforms was introduced and discussed. 


Chapter IV will present a MATLAB implementation model to identify and 
classify OFDM based wireless networking waveforms. This model is the foundation for 


all stimulations conducted throughout the thesis. 
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IV. BASEBAND SIMULATION RESULTS 


In this chapter, a MATLAB simulation is developed to implement the proposed 
waveform identification and classification scheme presented in Chapter III. The flow of 
the discussion will begin with a brief description of the simulation model. From this 
model, a MATLAB implementation is created, which will be the next topic discussed. 
Each section of the MATLAB implementation is covered in detail. Finally, the 
discussion will conclude with simulation results, interpretation of the results and 


conclusions. 


A. SIMULATION MODEL 


The proposed scheme outlined in Chapter III is implemented by MATLAB code 
as depicted in Figure 18. The simulation commences with the signal generator producing 
a single frame of either IEEE 802.11 or 802.16 compliant time samples. This frame will 
begin with the appropriate preamble followed by a predetermined number of OFDM 
symbols. The symbols will consist of user data and be constructed as described in 
Chapter II. Next the vector of time samples passes through an AWGN channel and 


proceeds to the preamble correlation section. 


The noisy time samples are subjected to a cross-correlation operation against 
standard compliant preamble samples for IEEE 802.11 and 802.16. The results of the 
cross-correlation operation are utilized to make a decision, based on predetermined 
criteria, as to whether an IEEE 802.11 or 802.16 waveform is identified. If the 
predetermined threshold is not met, the simulation ends. Once this decision is made, if an 
IEEE 802.11 signal is identified it is processed by the FAM to produce the SCD estimate 
and confirm the preamble correlation results. If an IEEE 802.16 signal is detected, the 
signal is processed by a test FAM operation to classify the CP length. The test FAM is a 
specially configured SCD estimator which will identify the CP length of the IEEE 802.16 
signal. Since the IEEE 802.11 standard only provisions for a quarter length CP, this 
process is not necessary for 802.11 signals. With out classification of the CP length in 


the case of IEEE 802.16, the FAM would not be able to extract the signal’s 
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cyclostationary signature. Once the cyclostationary signature has been produced, the 
signal is demodulated and compared to the transmitted data string. This, obviously, is to 
ensure the integrity of the transmission, identification and cyclostationary analysis 


processes. The next section examines in detail the main MATLAB program. 





802.11, 802.16 a 
Signal AWGN 
Generator Channel 





Figure 18. MATLAB Implementation Model. 


B. IMPLEMENTATION 


The MATLAB implementation of the model depicted in Fig. 18 is performed in a 
modular fashion. The main program provides a skeletal structure with individual 
operations being called as functions. All user-defined inputs reside within the main 
program and are imported to the function routines where applicable. Figure 19 is a flow 
chart representation of main program. The routine begins with the user selecting the 
desired waveform to be generated and the desired SNR in the channel. The appropriate 
waveform generator function is then called to produce one frame of either IEEE 802.11 
or 802.16 compliant OFDM symbols with preamble symbols at the beginning. Next the 
vector of waveform samples has AWGN added to it to simulate the channel. The 
variance of the AWGN is determined by the SNR selected by the user. At this point in 
the routine, the preamble of the noisy sample vector is cross-correlated with 320-sample 
long standard compliant IEEE 802.11 and 802.16 preamble sequences. By meeting or 
exceeding predetermined identification criterion, the preamble correlation identifies the 
waveform as being 802.11 or 802.16 compliant. If the identification criterion is not met, 


the routine terminates as no match is detected. 
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In the event an IEEE 802.11 waveform is identified, the program will call a FAM 
SCD estimation function, which will generate the SCD estimate of the waveform. The 
main program then analyzes the results of the SCD function to confirm the results of the 


preamble correlation operation. 


If an 802.16 waveform is identified, the routine will first classify the CP and then 


confirm the results of preamble correlation in a similar process as with 802.11 


waveforms. 
Select 
Waveform to be 
Generated 
Create 802.11 Call Wavela Create 802.16 
Waveform Generator Waveform 
BLAI CON! 

802.11 ‘ 802.18 
Preamble reamble Preamble 
Identified Correlate Identified 
Confirm Identify CP size 

via FAM Pilot No Match via FAM 
SC Signature End Routine Analyzer 
Demodulate Confirm 
802.11 via FAM Pilot 
Signal SC Signature 
Signal Classified aaa 
Signal 
Signal Classified 
Figure 19. Flow Chart of Main Program. 
iis Baseband Signal Generation 


In order to produce IEEE 802.11 and 802.16 compliant waveforms for simulation, 
two MATLAB programs were developed. Generation of each waveform is accomplished 


via WiFi_BasebandMod.m and WiMax _BasebandMod.m, two separate functions called 
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by the main program. The first function produces IEEE 802.11 compliant waveforms 
and the second produces IEEE 802.16 compliant waveforms. Each of these functions 
follows a similar pattern. Figure 20 illustrates the sequence of events that transpire when 
the signal generator function is called. First, the signal generator creates a symbol 
constellation for the appropriate baseband modulation scheme. The baseband modulation 
scheme used is determined by a logic sequence in the main program and is a function of 
the user-defined channel SNR. For an IEEE 802.11 waveform, a 0.8-ms down link frame 
is used. This equates to one preamble sequence, a signal symbol and 195 OFDM 
symbols being created, each 4 us in length. An IEEE 802.16 frame length is set at 10 ms 
but only the downlink portion of the frame is utilized. This equates to two preamble 


symbols and 67 control and data OFDM symbols. 
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Figure 20. Generic flow of Signal Generator. 


To form an OFDM symbol, first a column vector of baseband modulated I-Q data 
of size equal to the number of data subcarriers is formed. For IEEE 802.11 waveforms 
the size is 48 and for 802.16 it is 192. The pilot subcarriers are inserted with the 
corresponding standard defined pseudo-random sequence at the appropriate subcarrier 
frequency offset index locations. The guard subcarriers (nulls) are then included to 
complete the vector for the IFFT operation. After the IFFT operation, the CP is appended 
to the front of the IFFT output values, thus completing the formation of an OFDM 
symbol. 


Data manipulation, to include Forward Error Correction (FEC) coding, data 
scrambling, and interleaving, are intentionally omitted from both signal generators as 


they have no impact on the cyclostationarity of the waveform. 
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Each signal generator has the capability to generate BPSK, QPSK, QAM-16 and 
QAM-64 baseband modulated waveforms. Figure 21 illustrates the I-Q constellation of a 
transmitted QAM-64 modulated signal. Of note, the data point located at zero in-phase 
voltage and zero quadrature voltage represents the DC null. The I-Q values located at 
(-1,0) and (1,0) coordinates represent the pilot subcarrier BPSK modulated symbols. The 
remaining 64 I-Q values represent the data symbols carried by the data subcarriers. 
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Figure 21. I-Q Voltage Constellation of Quadrature Amplitude Modulation 64 
with DC/Guard Nulls and Pilot Symbol Locations. 


2. Preamble Cross-correlation Identification 


As described in Chapter II, IEEE 802.11 and 802.16 waveforms have their own 
distinct preamble sequences, which are transmitted at the beginning of each burst 
transmission. Because each waveform’s preamble is unique, it provides a reliable 
method of identifying the presence of a standard compliant waveform. There are two 
correlators in this portion of the main program that compare 320 samples of the received 


waveform with 320 samples of reference waveform. The first correlator’s reference 
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waveform consists of 320 samples of IEEE 802.11 standard compliant data. The second 
correlator’s reference waveform is composed of 320 samples of IEEE 802.16 standard 
compliant data which consists of the first OFDM preamble symbol with a 1/4 length CP. 
By using the 1/4 CP length option as the reference waveform, we ensure that the entire 
preamble sequence will be present in the received waveforms. If the received and 
reference waveforms match, the cross-correlation process will produce a clearly 
identifiable peak for /=0 in Equation 3.7. If no peak is present from either correlator, 
the received waveform is neither IEEE 802.11 nor 802.16 and the program will terminate. 


Identification criteria will be addressed in the results portion of the chapter. 


The next operation performed by the main program will depend on the waveform 
that is identified. If an 802.11 waveform is identified, the next step will be to confirm the 
identification by its cyclostationary signature. If the received waveform is identified as 


802.16, the CP will next be classified as described in the next section. 


3. 802.16 CP Classification 


This section explains the process in which the test FAM SCD function determines 
the CP size of the received signal. Although preamble correlation is an effective method 
of identifying the waveform received, an IEEE 802.16 signal can not be decoded without 
knowledge of the CP size. The length of the CP is determined when the network is setup. 
There are no transmissions that identify the CP length of an IEEE 802.16 waveform, so 
the passive listener needs a method to classify the CP length. IEEE 802.11 waveforms 
always use a CP that is 4 the length of the OFDM symbol. For this reason, IEEE 802.11 
identified waveforms are not subjected to the test FAM operation. A method of 
identifying the CP length is necessary in order to know how many samples to discard at 
the beginning of each OFDM symbol. In addition, when the CP is not identified and 
removed prior to cyclostationary analysis by the FAM estimator, the cyclostationary 


signature of the signal is significantly diminished. 


Figure 22 depicts a flow chart representation of the process that classifies the CP 
of IEEE 802.16 identified signals. After the preamble correlation operation identifies that 


an IEEE 802.16 waveform is present, the received signal vector is imported into the IEEE 
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802.16 autofamtest.m function. The IEEE 802.16 autofamtest.m function is a slightly 
modified version of the FAM SCD function estimator described in Chapter III. 
Parameters imported to the IEEE 802.16 autofamtest.m are the same as with the FAM 
function with the exception of the CP size. Since the CP length is unknown to the 


receiver, the IEEE 802.16 autofamtest.m is run for the four different possible CP sizes. 
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Figure 22. Program Flow to Classify IEEE 802.16 CP length. 
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The CP classification process begins with the IEEE 802.16 autofamtest.m 
processing the signal assuming it was generated with a 1/32 length CP. The received 
signal is imported into the IEEE 802.16 autofamtest.m function along with preselected 
values for the sampling frequency, frequency resolution, and cyclic frequency resolution, 
just to name a few. The first step performed by the function is to format the signal 
vector. This process begins by discarding the preamble samples of the received signal 
data vector. Within both IEEE 802.11 and 802.16 compliant waveforms, the preamble 
symbols are nonstandard OFDM symbols, as described in Chapter I. They would tend to 
diminish any cyclostationary pilot subcarrier signature that is being generated by the 
estimator. After discarding the preamble samples, the CP of each OFDM symbol is 
discarded by removing the first eight samples of every successive 264 sample blocks 
throughout the input data vector. This removes a 1/32 CP amount of samples from each 
symbol in preparation for FAM processing. Once this process has been performed on the 
entire vector, the data is processed by the test FAM estimator as mentioned in Chapter II 
with one notable exception. In order to develop the cyclic pattern of the pilot subcarriers, 
more than one OFDM symbol must be processed. By processing multiple symbols, the 
spectral lines produced by the pilot subcarriers’ cyclostationary signature will become 


more pronounced. 


The data subcarriers are modulated with random I-Q data, and thus, possess no 
deterministic reoccurring pattern. The pilot subcarriers are modulated with pseudo- 
random BPSK data and possess a small degree of repetition. In order to exploit and 
strengthen this marginal cyclic pattern, SCD results from multiple blocks of samples are 


averaged. 


The averaged SCD values from the FAM operation are displayed as shown in 
Figure 23, which depicts the principal cyclic frequencies along the x axis. When a robust 
pilot subcarrier cyclostationary signature is generated, identifiable spectral lines will be 
displayed along this axis. These spectral lines will be located at the respective pilot 


subcarrier position throughout the cyclic spectrum of the signal The shaded portion of 


the S@ ( Fi ) matrix in Figure 23 is the region of support for the SCD function of the 
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processed signal within the 257x1025 array (IEEE 802.16 SCD matrix dimensions). 


Values within this region represent SCD magnitude estimates. 
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257 Rows 
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Figure 23. Matrix Description of FAM SCD function estimator 
Output, S% (f). 


Since the row containing the principal cyclic frequency values has been identified, 
the individual cells containing the pilot subcarrier spectral lines now need to be located. 
As mentioned in Chapter I, IEEE 802.16 compliant waveforms employ a 256 point IFFT 
to create each OFDM symbol. With the principal cyclic frequency row having 1025 
cells, it was determined that floor(1025/4) or approximately 4 cells represented a 
subcarrier. This methodology is utilized to determine an approximate location of the 
eight pilot subcarrier spectral line locations. Next, the average of the magnitude of these 
eight spectral lines is determined. There is one cell in this row vector for each possible 
CP length. The CP length that generates the most robust cyclostationary signature will 


produce the largest average pilot subcarrier magnitude. 


Cyclostationary feature identification is now used to confirm the results generated 


by the preamble correlation process. 
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4. Cyclostationary Waveform Identification 


Although the received waveform was identified by preamble correlation, the next 
portion of the model will confirm that result through cyclostationary feature 
identification. In the case of a received IEEE 802.16 waveform, it has already been 
subjected to cyclostationary identification when the CP was classified. Now a technique 
similar in nature to the CP classification process will be utilized. The following 
discussion will be broken down into two sections. One section will cover the specifics of 


IEEE 802.11 cyclostationary identification and the second will describe 802.16 specifics. 


a. IEEE 802.11 Cyclostationary Feature Extraction 


This step of the implementation model is to confirm the preamble cross- 
correlation waveform identification results produced by the correlators using 
cyclostationary feature extraction. The process is quite similar to the CP identification 
discussed in the previous section. The function WiFi _autofam.m is called to produce the 
SCD function of the 802.11 waveform. The settings used to control the base line 
parameters of the FAM SCD estimator are imported when the function is called. Some of 
the parameters imported to the function include: signal vector famin, sampling 
frequency f., frequency resolution Af and cyclic frequency resolutionAa. As with the 
CP classification process, the signal vector data is first formatted by discarding the 
preamble and signal symbol samples at the beginning of the frame. Next the CP is 
removed from every OFDM symbol throughout the remainder of the frame. Then the 
formatted signal vector is processed by the FAM SCD estimator. The SCD function, 


Sé (f), is averaged over multiple OFDM symbols. 


When simulating the IEEE 802.11 waveform, an inter channel spacing 
factor of 20 MHz was used. Optimal values for Af and A@ where determined to be 
156,250 Hz and 78,125 Hz, which were 2 and % the inter subcarrier spacing distances, 


respectively. These parameters generated output array dimensions for S%(f ) of 


129x513. Row 65 of this array contained the principal cyclic frequency values, to 
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include four pilot subcarrier spectral lines. The magnitudes of the pilot subcarrier 
spectral lines are then averaged in the following fashion, 
1 
P tot os by Six (0) (4. 1) 
4 ie{pilot indices} 


where { pilot indices} = {-21, -7,7,21}. In addition to averaging the pilot subcarriers, 


the data subcarrier’s magnitudes are also averaged as shown 


1 
lone Shas 
48 ie{data indices} 


Si, (0) (4.2) 


where {data indices\ = {—26 :22,-20:-8,-6:-1,1:6,8:20,22:26}. The guard and DC 


subcarriers are omitted as their inclusion would have drastically altered the data 


subcarrier average. 


Now that the average peak magnitudes for the pilot and data subcarrier 
cyclic spectral lines have been determined, they are compared. The metric used for 
classification is threshold y given by 


Pe 

pilot 

= er 4.3 
i (4.3) 


data 
For classification implementation a 10% margin was chosen leading to the following 
classification criteria 
Hoo ysl 
ee (4.4) 
gt vel A 


where H, is a positive classification of an IEEE 802.11 waveform and H, is a null 


hypothesis. 

The next section will discuss the IEEE 802.16 specifics of cyclostationary 
feature extraction. 

b. IEEE 802.16 Cyclostationary Feature Extraction 


The process of cyclostationary feature extraction employed on 802.16 
compliant waveforms is the same as the 802.11 process with one exception; different 


settings are used to control the base line parameters of the FAM SCD estimator, IEEE 


Sil 


802.16 autofam.m. IEEE 802.16 simulations utilized a transmission bandwidth of 3.5 
MHz. This bandwidth required a sampling frequency of 4.0 MHz, which is determined 
in the following fashion, 
8 W 
f= ea 8000 (4.5) 
where W is equal to the bandwidth of the signal. The optimal frequency and cyclic 
frequency resolution were determined to be the inter subcarrier spacing of 15,625 Hz for 


Af and half this value, 7813 Hz, forA@. As with the 802.11 process, the row containing 


the principal cyclic data is identified and the location of the eight pilot subcarrier spectral 
lines identified. These magnitudes were averaged as follows 
Pitot oe > vee (0) > (4.6) 
8 ie{pilot indices} 


where {data indices} = {—88,—63, -38,—13, 13,38, 63,88}. The data subcarriers where 





averaged in a similar fashion as shown 


l . 
P, =——= si(0), 4.7 
e 192 os : ( ) ) 





where {data indicies\ = { 100 : -89, -87 : -64, -62 : -39, -37: -14,-12: -1,1:12,14:37,39:62 


64:87,89:100}. These values are subjected to the same classification criteria discussed 
in the previous section. 


A brief discussion of the receiver section of the implementation model 


will now be conducted. 


5. Baseband receiver 


This step of the implementation model is to ensure the data transmitted is 
recoverable. Although this step is not vital to the identification and classification 
process, it does lend credibility to the signal generator and allows for the observation of 
the channel effects on the signal. Both IEEE 802.11 and 802.16 data recovery operations 


are similar and will be briefly reviewed. 
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The receiver operation is relatively simple compared to the transmitter operation. 
The first step the receiver performs is to discard the preamble. Next, in a symbol by 
symbol fashion, the CP is removed. At this point, it becomes apparent why proper CP 
classification is critical for an 802.16 waveform. Without classification of the CP, the 
receiver would not know the proper number of samples to discard for each symbol CP 
and the data represented by the OFDM symbol samples would become unrecoverable. 
The I-Q data is now processed by an FFT operation in preparation for baseband 
demodulation. Next the guard and pilot subcarrier data is discarded. Once the baseband 
demodulation is performed, the recovered binary sequence is compared to the original to 
determine bit errors. Figure 24 is an illustration of a received QAM 64 constellation with 
the channel introducing a SNR of 20 dB. A bit error rate of 0.004618 was determined 
after data recovery and comparison in this case. For comparison, the theoretical 


probability of bit error is calculated from [14], 





Pi, 2 


1 
eae 
[ VM 3log, VM 2E, 
QO (4.8) 
log, (VM } Cat) ON, 

where M is equal to 64 and E£;/N, is equal to 20 dB. The result of the calculation 
yields3.7x10~, which is reasonably close considering that Equation 4.8 is an 
approximation and the sample size of the simulation was 156,672 bits long, not an 


adequately large sample size. 
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Figure 24. Received I-Q Voltage Constellation of Quadrature Amplitude 
Modulation 64 with Channel Conditions Simulating a SNR of 20 dB. 


The next section with cover the results of the implementation model, as well as a 


few simulations to help identify some limitations of the model that we need to be aware 


of. 


C. RESULTS 


The discussion now shifts to the results of the implementation model. The first 
results addressed will be the preamble correlation findings, followed by CP classification 


and finally the cyclostationary feature extraction. 


1. Preamble Cross-Correlation 


Figures 12 and 13 from Section A of Chapter III display the output of each 
correlator with both IEEE 802.11 and 802.16 waveforms passed through them at an SNR 


of zero dB. This process was extended to multiple SNRs in order to identify when the 
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results would become overwhelmed by the power noise and to identify a point at which a 
decision threshold could be set. Figure 25 illustrates the results of a simulation run to 
determine the minimum SNR that still produced reliable results from each correlator 
when an IEEE 802.11 preamble sequence is received. In addition to running this 
simulation at 21 different SNRs, all SNR iterations were repeated 150 times in order to 
average out the effects of noise. Two input waveforms were averaged by each correlator 
for this simulation, the received 802.11 waveform and a white Gaussian noise only 
waveform. Although the 802.11 correlator produces the largest average peak at all SNRs 
tested, a significant separation between this curve and the others occurs at —14 dB, where 
the matched correlator output approaches a steady value of 0.01256 while the 
uncorrelated output continues to decrease towards 0.00145. The correlation with white 


Gaussian noise input waveform approaches zero in both correlator outputs. This agrees 


with the results of Equation 3.7 where s, (t) is noise only. 
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Figure 25. Preamble Correlation Results of an IEEE 802.11 Received 


Waveform versus SNR Averaged Over 150 Runs. 
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Figure 26 illustrates the results of the same simulation when an 802.16 waveform 


or noise only waveform is received. As with the 802.11 results, a significant separation 


between the 802.16 correlation results and the remaining three results occurs at —14 dB. 


The 802.16 results approach a value of 0.01011, with the matched correlator output, 


while the unmatched results decrease to 0.001931. Again, the noise only waveform 


results decrease toward zero correlation as SNR increases. 
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Figure 26. Preamble Correlation Results of an IEEE 802.16 Received 
Waveform versus SNR Averaged Over 150 Runs. 


The purpose of this simulation was to demonstrate that preamble cross-correlation 


identification is effective at determining which waveform is received at an SNR range 


well below the receiver’s capability to recover the data. Based on visual observation of 


the simulation results displayed in Figures 25 and 26, without a detailed application of 


detection theory, a threshold of 0.009 was chosen as the decision metric for waveform 
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identification. By choosing 0.009 as the threshold, reliable identification decisions would 


be a reasonable expectation down to an SNR of —14 dB in both cases. 


Next, the discussion moves on to the results of the CP classification portion of the 


simulation. 


2. IEEE 802.16 CP Classification 


In order to validate the CP classification operation discussed in Section B-3, a 
simulation was performed to determine the minimum number of OFDM symbols that 
were required to be processed for a correct classification. There are four possible 
outcomes of this routine, one for each CP length that is possible: 1/4, 1/8, 1/16, 1/32. The 
test FAM SCD estimator tested a 1/4 length CP waveform as described in Section B-3 
against all four possible outcomes 200 times. The results were obtained for an integer 
number of OFDM symbols from 3 to 24 symbols. To measure the effectiveness of the 
classification, the metric used is CP classification percentage: number of successful tests 
out of the total number of trials (200 in this case). Figure 27 illustrates the results of this 
simulation. Notice, by 20 OFDM symbols, the test FAM will classify the correct CP 100 


percent of the time. 
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Figure 27. CP Classification Percentage versus the Number of OFDM 
Symbols Processed by the Test FAM SCD Function Estimator. 
Averaged Over 200 Runs. 
3. Cyclostationary Feature Extraction 


The final step of the implementation model for both waveforms is to confirm the 
identification and CP classification operations by cyclostationary feature extraction. The 
cyclostationary signatures of each waveform are presented as well as the results of two 
simulations conducted to determine minimum threshold requirements. The first 
simulation is designed to identify the minimum SNR where the cyclostationary feature 
extractor still produces reliable results. The second simulation is designed to determine 


the minimum number of OFDM symbols that needs to be processed by the FAM SCD 


estimator to produce reliable results. 


a. IEEE 802.11 


Figures 28-30 illustrate the cyclic signature exploited by cyclostationary 
feature extraction. These figures are generated using the MATLAB function 
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WiFi autofam.m. Figure 28 is a surface plot of the positive cyclic frequencies of 


Se ( f ). The left portion of the graph represents the PSD of the waveform with the 


positive cyclic frequency realizations extending to the right. The pilot subcarrier peaks 
can be seen protruding well above the data subcarrier floor in arrays of three then two and 
finally one. To ensure the pilot subcarrier peaks are clearly identifiable, the results 


averaged over 120 OFDM symbols. 
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Figure 28. Surface Plot Representation of an IEEE 802.11 Waveform’s 
S°(f), Averaged Over 120 OFDM Symbols. 


Figure 29 consists of the profile plots that display SCD versus the 
principal cyclic frequency axis and PSD versus the principal frequency axis. The pilot 
subcarriers spectral lines are not identifiable in Figure 29(a) and clearly discernable in 


Figure 29(b). 
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Figure 29. Profile Plots: Magnitude of S% (f) versus a) Cyclic Principal 


Frequency and b) Frequency for IEEE 802.11 Waveform. Results 
Averaged Over 120 OFDM Symbols. 


Figure 30 provides contour plot of S®% ( f ) at a level of 0.5. Unlike 


Figure 28, the contour plot includes the negative cyclic frequencies. The PSD is clearly 
identifiable, located at the zero cyclic frequency range with the DC null located at the 
center of the graph. As the cyclic frequency components increase in magnitude, the pilot 


subcarrier locations appear in a symmetric fashion. 
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Figure 30. Contour Plot of an IEEE 802.11 Waveform’s S% ( f) Magnitude. . 
Results Averaged Over 120 OFDM Symbols. 


Figures 31 and 32 were generated from simulations that were intended to 
determine the effects of the SNR and number of OFDM symbols on the subcarrier SCD 
estimates. Figure 31 shows the effect of SNR on the peak pilot and data subcarrier 
magnitudes. The results are obtained by averaging results from 150 simulation runs. At 
roughly -15 dB, the pilot subcarrier SCD peak value clearly began to emerge from the 


data subcarrier values. 


Figure 32 shows the effect of the number of OFDM symbols used on the 
SCD estimates of pilot subcarriers and data subcarriers. The results are intended to 
determine how many OFDM symbols are required to produce a clearly identifiable 


separation between the pilot and data subcarrier SCD peaks. 
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Figure 31. IEEE 802.11 Subcarrier SCD Values versus SNR. 
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Figure 32. IEEE 802.11 Averaged Subcarrier SCD Values versus Number of 


OFDM Symbols Processed by the FAM SCD Function Estimator. 
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b. IEEE 802.16 


Plots similar to those obtained for 802.11 waveforms were generated for 
an 802.16 waveform: a surface plot, a profile plot and a contour plot. These plots were 
generated from a waveform possessing a 1/32 length CP with results averaged over 120 
OFDM symbols using the MATLAB function IEEE 802.16 autofam.m. Figure 33 
depicts seven rows of pilot subcarrier cyclic spectral lines. Of these, only four possess 
pilot subcarrier spectral lines that reside along the principal cyclic axis (f= 0). The peak 
magnitude of the averaged SCD values are normalized to one. Further detail of the 


principal axis is provided in Figure 34. 





Sf) 





0 
2 25 3 BS 

‘ 

a (Hz) x10 


Figure 33. Surface Plot Representation of an IEEE 802.16 Waveform’s 
S“(f). Results Averaged Over 120 Symbols. 
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Figure 34 shows the profile plots of the SCD magnitude as a function of 
the cyclic frequency @ and frequency f. Figure 34(a) clearly shows all the pilot 


subcarrier peaks along the principle cyclic axis. 
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Figure 34. Profile Plots: Magnitude of S% (f) versus a) Cyclic Principal 


Frequency and b) Frequency for IEEE 802.16 Waveform. Results 
Averaged Over 120 Symbols. 


Figure 35 provides the top down perspective of the pilot subcarrier cyclic 
peak locations. The contour level in Figure 34 was set at 0.54 to ensure all the pilot 
subcarrier cyclic spectral lines were represented, while excluding spurious responses due 


to random data on data subcarriers. 
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Figure 35. Contour Plot of an IEEE 802.16 Waveform’s S@ ( .s ) Magnitude. 
Results Averaged Over 120 Symbols. 


Results of Figures 36 and 37 are similar to those in Figure 31 and 32, 
respectively. Figure 36 demonstrated a clearly identifiable separation at -10 dB between 
the pilot and data subcarriers. The results displayed in Figure 37 illustrate how the 
cyclostationary signature of the pilot subcarriers is not identifiable until more than three 
symbols have been averaged. In order to satisfy the above mentioned classification 
criteria, roughly 10 OFDM symbols need to be averaged. Beyond 20 symbols, not much 


improvement was obtained. 
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Averaged Subcarrier SCD Magnitude 
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Figure 36. IEEE 802.16 Averaged Subcarrier Peak Values versus SNR. 
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IEEE 802.16 Averaged Subcarrier Peak Values versus Number of 
OFDM Symbols Processed by the FAM SCD Function Estimator. 
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In summary, the surface, profile and contour plots of the previous section 
are clear illustrations of the effectiveness of the cyclostationary feature extraction 
techniques in regard to identifying and classifying OFDM based waveforms. Although 
the plots were generated by averaging 120 OFDM symbols, the implementation model 
coding is able to reliably identify the pilot subcarrier spectral peaks in as few as 20 
OFDM symbols. In both cases, after 20 OFDM symbols are averaged, the pilot 


subcarriers always possessed a higher magnitude than the data carriers. 
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Vv. CONCLUSIONS 


This thesis discussed the fundamentals of OFDM and how OFDM modulation is 
implemented in the IEEE 802.11 and 802.16 standards. Cyclostationarity was then 
introduced to provide the reader with a basic understanding required for the following 
chapters. Next, a scheme for identifying and classifying 802.11 and 802.16 waveforms 
was proposed and discussed. Three possible cyclostationary feature extraction 
approaches were discussed and vetted for implementation. The most promising method, 
pilot subcarrier cyclostationary feature extraction, demonstrated the most potential and 
was chosen to classify the CP size of 802.16 waveforms. Additionally, pilot subcarrier 
cyclostationary feature extraction was employed to confirm the results of preamble cross- 
correlation identification process. Finally, a MATLAB simulation model was developed 


and the simulation results presented and discussed. 


A. SIGNIFICANT RESULTS AND CONTRIBUTIONS 


There were four major contributions made by this thesis. First, the work 
established an approach to identify a received waveform with no a priori knowledge of its 
origin. After much research and experimentation, it was determined that a preamble 
cross-correlation operation would provide the most reliable results while minimizing 


computational complexity of the simulations. 


Next, a method of classifying the cyclic prefix length of IEEE 802.16 waveforms 
was devised. Since the CP option is established when the network is first setup, there is 
no way for a passive listener to identify this value by decoding captured control 
transmissions. By employing cyclostationary feature extraction, we were able to 


determine the CP length in a rapid fashion. 


After identifying three methods of implementing cyclostationary feature 
extraction from [1], [3] and [4], we needed to determine which method or methods would 
be most compatible with the stated objectives of this work. Although embedded 
cyclostationary signature extraction is a promising technique, pilot subcarrier 
cyclostationary feature extraction proved to be the most applicable method. Unlike the 
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waveforms analyzed in [1], IEEE 802.11 and 802.16 waveforms did not amplify their 
pilot subcarrier transmissions so an averaging process had to be devised and 


implemented. 


Finally, an overall scheme was developed to integrate the previous three processes 


into an organized framework for implementation. 


B. FUTURE WORK 


Initial identification of a received waveform was accomplished through a 
preamble correlation process. Although this method proved to be reliable, a rigorous 
application of detection theory was not applied due to time constraints. In order to 
generalize the results of the preamble correlation process, rigorous signal detection 


methods could be applied to this portion of the work to lend it greater credibility. 


In this thesis, cyclostationary feature extraction was used to identify and classify 
OFDM based waveforms. The construction of the OFDM symbols in the two IEEE 
standards referenced in this thesis is a static process in that the format of the symbols 
remained the same. Other than in the preamble symbols, the null, pilot and data 
subcarrier locations remained the same. This is not the case with the IEEE 802.16e 
standard, which employs Orthogonal Frequency Division Multiple Access (OFDMA). 
OFDMA combines the baseband modulation operation and the multiple access process 
into one operation. | Users are assigned a group of subcarriers, i.e., the number and 
location of data and pilot subcarriers will now vary with the number of users. So not only 
is the sequence that modulates the pilot subcarriers pseudo random, but with 802.16e, the 
location of many of the pilot subcarriers is also pseudo random, with respect to time. In 
order to obtain the pilot subcarrier cyclostationary signature, the Down Link MAP (DL- 
MAP) messages must be decoded to determine the number and location of both the fixed 


and non-fixed pilot subcarriers. 
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APPENDIX 


This Appendix contains the MATLAB code used in the implementation model. 


The main program is listed first with the functions called by this program listed in order 


of their utilization. 
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Select Waveform to be Created 


Ty 


WiFi 


1, 


WiMax=2 
Waveform_T 


Enter WiMax Or WiFi to create waveform 
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ype=WiMax; 


Signal to Noise Ratio of AWGN channel 
Logic based on SNR to determine Baseband Mod 
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SNR=18.1 
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m=1; 
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&& 


((SNR > 8) 
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Calls WiFi tx funct 
feeds required data 
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Preamble _corr2 


Funt with sample Preample 


Vect with first 320 samples 


,CPsym, 'biased') 
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Cross Correlation") 


Preamble_corr 
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xcorr ( 


Comp 


-015]) 
plot (xaxis,WiMax_Norm) 


802.16 Preamble/Recived Signal 


bel('Sample Shift') 
2015] > 


plot (xaxis,WiFi_Norm) 
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corr2,lags2] 
Preamble_corr 
abs ( 
linspace (-319, 319, 639) 


corr, lags] 
[WiFi _Peak,WiFi_Peak_Pos] 


320,1) 
abs ( 
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subplot (2,1,2), 


title('b) 
-330 330 0 


y (1 
hold off 


{Preamble 





WiMax _Norm 
bel ('Normalized Preamble Correlation") 


bel ('Normalized Preamble Correlation") 


bel('Samples Shift') 








[WiMax_Peak, WiMAx Peak_Pos] 


[WiMax _Pre,WiFi_Pre] 
WiFi _Norm 
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subplot (2,1,1), 
yla 


CPsym 
xaxis 
hold on 
figure (7) 
axis ( 
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-009 & WiMax Peak > 


if WiFi_Peak > 


num2str('System Error, 
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Both Waveforms Identified") 


else if WiFi_Peak < 0.009 & WiMax Peak < 0.009 


num2str('Unknown Waveform") 
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else if WiFi_Peak > 0.009 


Estimate 





Sampling Freq used by FAM 
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20000000; 


fs= 


Freq Resolution of FAM 


ceil (20*10%6/128) 
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ceil (20*10%6/256); % Cyclic Freq Resolution of FAM 


dalpha= 


ss 


Eh CP 


ded into FAM Estimator 


10a 





1 wi 
Sets increment size of OFDM symbols for FAM 


Size of OFDM symbo 
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Transposes tx signal+AWGN vector 


WiFi_autofam 


% 


=y'; 


famin 


(famin,fs,df,dalpha,ss, index,v, 


[Sxp,fo,Np, centrow, centcol] 


shift) 


sum([Sxp(centrow, 89),Sxp(centrow,201),... 
Sxp (centrow, 313),Sxp(centrow, 425) ])/4 


PilotSC_ave 


200) pews 
312),Sxp(centrow, 321 


sum([Sxp (centrow, 49:88) ,Sxp(centrow, 97: 


Sxp (centrow, 209 


DataSC_ave 


424) pea 


248),Sxp(centrow, 265 


464) ])/376 


if PilotSC_ave > DataSC_ave*l1.1 


Sxp (cent row, 433 


num2str('IEEE 802.11 Waveform Classified') 


num2str('IEEE 802.11 Waveform Not Classified") 
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Calls rec funct 
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(L,N,m,Ndata,CP,y,bk, Frame) 





WiFi _BasebandDemod 


[ck] = 


feeds required 


Adjusts length of input data for 
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length (ck) ) 


bkl=bk (1,1 


Compares rec data string to trans 
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Error Rate 





Calculates Bit 


max(cumsum(err))/(m*Ndata*N) &% 


R= 


= 
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Input settings for WiMax FAM Spectral Correlation Density Function’ 
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else if WiMax _Peak > 0.009 


Sampling Freq used by FAM Estimate 


Freq Resolution of FAM 


% 


4000000; 


fs= 


% 


ceil (15625); 


df= 


73 


dalpha=ceil (15625/2) ; 
ss=CP+L; 

FAM 

index=30; 

analysis 68 


ae 


Cyclic Freq Res of FAM 
Size of OFDM Sym w/CP. Used to prep data for 


oe 


ol? 


Sets increment size of OFDM sym for FAM 


shift=512; & Num of data samples loaded into FAM each 
iter 
v=[.55]; & Height of contour plot plane 
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famin=y'; 
for num=1:4 








CP1=2% (num+2) ; % Sets CP length for famtest 

[Sxp, fo, Np, centrow, centcol]=WiMax_autofamtest (famin,fs,df,... 
dalpha,ss,index,CP1,shift) ; % Generates SCD for compison 

peak1 (1,1)=(Sxp(centrow,163)); 

peak1 (1,2)=(Sxp (centrow, 263) ); 

peak1 (1,3) =(Sxp (centrow, 363)); 

peak1 (1, 4)=(Sxp (centrow, 463) ); 

peak1 (1,5)=(Sxp (centrow, 563)); 

peak1 (1, 6) =(Sxp (centrow, 663)); 

peak1 (1, 7)=(Sxp (centrow, 763)); 

peak1 (1, 8) =(Sxp (centrow, 863)); 

peak (1, num) =mean (peak1); % Determines mean of peakl 

end % vector 


[Sigtype, location] =max (peak) ; 
if location== 
CP2=256/32 
Signal_Type=num2str('256 FFT with a 1/32 Cyclic Prefix') 
else if location== 
CP2=256/16 
Signal_Type=num2str('256 FFT with a 1/16 Cyclic Prefix') 
else if location== 
CP2=256/8 
Signal_Type=num2str('256 FFT with a 1/8 Cyclic Prefix') 
else if location== 
CP2=256/4 
Signal_Type=num2str('256 FFT with a 1/4 Cyclic Prefix') 














end 
end 
end 
end 


[Sxp, fo, Np, centrow, centcol]=WiMax_autofam(famin, fs,df,dalpha,ss, index,v 


peee 


shift; CPZ); % Creates SCD and graphs for the CP signal 
PilotSC_ave=sum([Sxp(centrow, 163),Sxp(centrow,263),... 


Sxp (centrow, 363),Sxp(centrow, 463),Sxp(centrow, 563), 
Sxp (centrow, 663),Sxp(centrow, 763),Sxp(centrow, 863) ])/8 
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rror Rate end 
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hm 





and feeds required data 
Adjusts length of input data for compare 


% 


Compares rec data string to trans 
Calculates Bit 


% 


fo) 
oO 


(L,N,m,Ndata,bk,CP,G,Frame) ; 


(1.1) 
num2str('IEEE 802.16 Waveform Classified') 


° 


(L,N,m,Ndata,CP,y,bk,Frame);sCalls rec funct 


E 802.16 Waveform Not Classified') 


Dp 
G 


Bt 





Written by Steven Schnur 


Baseband Data Setup and Modulator 


WiMax Receiver Section 
Generates 64QAM Modulation Object 


num2str('IE 


WiFi_BasebandMod 


sum([Sxp (centrow,117:160),Sxp(centrow,165:260),... 


Sxp (centrow, 265:360),Sxp(centrow, 365:460),Sxp(centrow, 465:512),... 


Signal_Type 


[x] 





WiMax_BasebandDemod 


Sxp (centrow, 517:560),Sxp(centrow, 565:660),Sxp(centrow, 665:760),... 


Sxp (centrow, 765:860),Sxp(centrow, 865:913)])/761 


if PilotSC_ave > DataSC_ave* 
error ('Wrong number of arguments") 


Signal_Type 
bk (1,1:length (ck) ); 
max (cumsum(err) ) / (m*Ndata*N) 


ER 


5 


hm 





DataSC_ave 
else 

end 

[ck]= 

bk1 

B 

function 

if nargin~=8 
end 
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Normalization Value of 64QAM 
(=7=19) SG. (TFT DG eas 
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IQ data 


Modulates m sized sym with obj: 
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modulate (object,mk') 


Xp= 


Stacks IQ data in a column vector 
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Ndata loop 
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End of w 


% 


end 


oe 
oe 
oe 
oe 
oe 
o\? 
oe 
oe 
oe 
oe 
oe 
a? 
ole 
oe 
le 
oe 
oe 
oe 
ole 
oe 
ale 
oe 
le 
oe 
oe 
oe 
le 
oe 
ale 
a? 
le 
oe 
oe 
oe 
le 
oe 
ale 
oe 
le 
oe 
ale 
oe 
le 
oe 
ol? 
oe 
ol? 
oe 
oe 
oe 
oe 
oe 
ale 
oe 
le 
oe 
oe 
oe 
ole 
oe 
ale 
oe 
le 
oe 
oe 
oe 
le 
oe 
ole 
a? 
ol? 


Ale ol? 


Pilot SC and Data in 


DC=0, 
appropriate SC locations 
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DC Null 
Data 


% 


0; 


X(1,1) 
X (2 


T 


30,1) 


7,1)=X1(25: 


Positive Pilot SC 


Data 


% 


¥ 


(1,2) *G*1 


Pilot_SC 


X(8,1)= 
X(9 


243, 1) ¥ 


21,1) =X1 (31 


Positive Pilot SC 


Data 


% 


=Pilot_SC (1,2) *G*-1,; 


X (22,1) 
X23 


748),.1:) 3 


27,1) =X1 (44 
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X(28:32,1)=0; % Lower Guard Band 
X(33:38,1)=0; % Upper Guard Band 
X(39:43,1)=X1(1:5,1); % Data 
X(44,1)=Pilot_SC (1,2) *G*1; % Negitive Pilot SC 
X(45:57,1)=X1(6:18,1); % Data 
X(58,1)=Pilot_SC (1,2) *G*1; % Negitive Pilot SC 
X(59:64,1)=X1(19:24,1); % Data 
Pilot_SC=circshift (Pilot_SC, [0,-1]);% Pseudo-randomize pilot seq IAW 
Xplot ((k-1) *L+1:k*L,1)=X; % WiFi standard 

end %& End of else if index(1,k)==1 loop 





2.9. 9. o- 























xt=ifft (X,L); % L size IFFT to convert I-Q data to time samples 
g(1:CP,1)=xt ((L-CP+1):L,1); % Loads CP worth of Data end of OFDM sym 
x1(1:CP,1)=g; % Loads CP to front of vector 

x1 (CP+1:L+CP,1)=xt; % Loads rest of OFDM Symbol into vector 
end % End line 97 loop 

end & End line 95 loop 

end % End line 93 loop 

end & End line 90 loop 

x (((k-1) * (L+CP) +1) :k* (L+CP),1)=x1; % Fills tx vector w/CP appended 
end % End for k=1:N loop 

figure (1) 

plot (xX, '*"); & Plot Transmitted Constellation 
title('Transmitted I-Q Data") 





xlabel('Transmitted In Phase Data') 
ylabel('Transmitted Quadrature Data') 
axis.([=1.5 1.5 -1.5 1.5]); 

grid on 
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COSC SCSC SSC SCCCSCCESCSCCC CCE CSCC CCE CSCC ECE CECCC CCE CCST ECE CECE CEEEC CCE SEEEC CS SSSES 
° . 4 ° 
6 WiFi Short Preamble Time Samples % 
[oa oan © Sa oa © al © Sa © al © Sa © Sam © al © Wa © Sas © al © Sa © al © Sa © Sas © al © Sa © Sal © al © al © Sal © al © Sal © al © al © Sa © al © Sa © Sas © al © a © Sal © al © Sa © Sal © a © Dal © al © al © Dal © al © Sa © Sal © al © Ua © al © al © a © al © Sal oa © a © al oa © eal © al © aL © alo Sal © ao al oa oa oe oo eo eo eaeD 
CCCCCCC CC CSCC CCC CSCC SC SSC SSC SC SSC SSCS SC SSCS SSCS SSC SSCS CSSCSCSSCSSCSSCSSCSSCS CSCS CSCS COS 


function [short1l, short2]=WiFi_short_preamble; 

















short1=[0.023+0.0233;-0.132+0.0023;-0.013-0.0793;0.143-0.013j;... 
0.092+0.0003;0.143-0.0135;-0.013-0.0793;-0.132+0.00235;... 
0.046+0.0463;0.002-0.1325;-0.079-0.0133;-0.013+0.143;... 
0.000+0.09235;-0.013+0.1435;-0.079-0.01335;0.002-0.132j;... 
0.046+0.0463;-0.132+0.00235;-0.013-0.0793;0.143-0.013j;... 
0.092+0.0003;0.143-0.0133;-0.013-0.0793;-0.132+0.00235;... 
0.046+0.0463;0.002-0.1325;-0.079-0.0133;-0.013+0.143j,;... 
0.000+0.0923;-0.013+0.1435;-0.079-0.0133;0.002-0.132)5; 
0.046+0.0463;-0.132+0.0023;-0.013-0.0793;0.143-0.013j;... 
0.092+0.0003;0.143-0.0133;-0.013-0.0793;-0.132+0.00235;... 
0.046+0.0463;0.002-0.13235;-0.079-0.0133;-0.013+0.143j;... 
0.000+0.09235;-0.013+0.1435;-0.079-0.01335;0.002-0.132j;... 
0.046+0.0463;-0.132+0.0023;-0.013-0.0793;0.143-0.013j;... 
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.092+0.0004;0.143-0.0134;-0.013-0.0794;-0.132+0.0024;... 
.046+0.0464; 0.002-0.1323;-0.079-0.0134;-0.013+0.1434;... 
.000+0.0924;0.013+0.1434;-0.079-0.0134;0.002-0.1323;... 
j;-0.132+0.00245;-0.013-0.0793;0.143-0.0134;... 
.092+0.0004;0.143-0.0134;-0.013-0.0794;-0.132+0.0024;... 
.046+0.0464;0.002-0.1323;-0.079-0.0134;-0.013+0.1433 
.000+0.0924;-0.013+0.1434;-0.079-0.0134;0.002-0.1323 
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short 2=[0.046+0.0463;-0.132+0.0024;-0.013-0.0793;0.143-0.0134;... 
.092+0.0004;0.143-0.0133;-0.013-0.0794;-0.132+0.0024;... 
.046+0.0464;0.002-0.1323;-0.079-0.0133 
.000+0.0924;-0.013+0.1434;-0.079-0.013 
+0.0465;-0.132+0.0024;-0.013-0.079 
+0.0003;0.143-0.0134;-0.013-0.0793 
+0.0465;0.002-0.1324;-0.079-0.0133 
+0.09245;-0.013+0.1434;-0.079-0.013 
+0.0463;-0.132+0.0024;-0.013-0.079 
+0.0003;0.143-0.0134;-0.013-0.0793 
+0.0465;0.002-0.1324;-0.079-0.0133 
+0.09245;-0.0134+0.1434;-0.079-0.013 
+0.0463;-0.132+0.0024;-0.013-0.079 
+0.0003;0.143-0.0134;-0.013-0.0793 
+0.0465;0.002-0.1324;-0.079-0.0133 
+0.0925;-0.013+0.1434;-0.079-0.013 
+0.0463;-0.132+0.0024;-0.013+0.079 
+0.0003;0.143-0.0134;-0.013-0.0793;-0. or 
+0.0465;0.002-0.1324;-0.079-0.0134;-0.013+0.1434;... 
.000+0.0924;-0.013+0.1434;-0.079-0.0134;0.002-0.1323]; 
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function [longl, long2]=WiFi_long_preamble; 
long1=[-0.078+0.0003;0.012-0.09835;0.092-0.1063;-0.092-0.115j;... 
—0.003+0.0545;0.075+0.0745;-0.127+0.02135;-0.122+0.01735;... 
—0.0354+0.1515;-0.056+0.0225;-0.060-0.08135;0.070-0.014j;... 









































0.082-0.09243;-0.131-0.0654;-0.057-0.0394;0.037-0.0984;... 
0.062+0.0623;0.119+0.00445;-0.022-0.1615;0.059+0.0154;... 
0.024+0.05945;-0.1374+0.0475;0.001+0.1154;0.053-0.0044;... 
0.098+0.0263;-0.038+0.1064;-0.115+0.0554;0.060+0.0884;... 
0.021-0.0284;0.097-0.0834;0.040+0.1114;-0.005+0.1204;... 
0.156+0.0003;-0.005-0.1204;0.040-0.1115;0.097+0.0834;... 
0.021+0.0285;0.060-0.0884;-0.115-0.0554;-0.038-0.1064;... 
0.098-0.0265;0.053+0.00445;0.001-0.1154;-0.137+0.0474;... 
0.024-0.0594;0.059-0.0154;-0.022+0.1613;0.119-0.0044;... 
0.062-0.0624;0.037+0.0984;-0.057+0.03945;-0.131+0.065;... 
0.082+0.0924;0.070+0.0144;-0.060+0.08135;-0.056-0.0224;... 
-~0.035-0.1514;-0.122-0.0173;-0.127-0.0213;0.075-0.0744;... 
~0.003+0.0544;-0.0924+0.1154;0.092+0.1064;0.012+0.0985;... 








-0.156+0.0004;0.012-0.0984;0.092-0.1064;-0.092-0.1153;... 
-0.003-0.0545;0.075+0.0745;-0.127+0.02145;-0.122+0.0174;... 
-0.035+0.1514;-0.056+0.0224;-0.060-0.0814;0.070-0.0144;... 
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0.082-0.0924;-0.131-0.0654;-0.057-0.0394;0.037-0.0983]; 
long2=[0.062+0.0623;0.119+0.0043;-0. 
.024+0.0594;-0.137+0.04745;0.0014+0.1 
.098+0.0264;-0.038+0.1064;-0.115+0.0554;0.060+0.088 
.021-0.0284;0.097-0.0833;0.040+0.1114;-0.005+0.1203 
.156+0.0003;-0.005-0.1204;0.040-0.1115;0.097+0.0833 
.021+0.0284;0.060-0.0883;-0.115-0.0554;-0.038-0.106 
.098-0.0264;0.053+0.0044;0.001-0.1154;-0.137-0.0473;... 
.024-0.0594;0.059-0.0154;-0.022+0.1614;0.119-0.00435;... 

.062-0.0624;0.037+0.09845;-0.057+0.0394;-0.131+0.0654;... 
.082+0.0924;0.070+0.01445;-0.060+0.0814;-0.056-0.0224;... 
-0.035-0.1514;-0.122-0.0173;-0.127-0.0213;0.075-0.074 
-~0.003+0.0544;-0.0924+0.11545;0.092+0.1064;0.012+0.0985;... 
-~0.156+0.0004;0.012-0.0985;0.092-0.1064;-0.092-0.1154;... 
-0.003-0.0544;0.075+0.0745;-0.127+0.0214;-0.122+0.017 
-0.035+0.1514;-0.056+0.0224;-0.060-0.0814;0.070-0.014 


022-0.16145;0.059+0.0153;... 
154;0.053-0.0044;... 
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50.119+0.0044; 
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-~0.022-0.1613; 


4;0.001+0.1154; 
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70.037-0.0983 


0.059+0.0154;... 


0.053-0.0043;... 
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.021-0.02835;0.097-0.0835;0.040+0.1115;-0.005+0.12035]; 
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function [WiMax_Pre,WiFi_Pre]=Comp_Pre; 

WiMax _Pre=1*... 

[-0.0,0.0997 — 0.0305i1,-0.0055 + 0.01401,0.0949 + 0.1927i,... 

0.0735 + 0.0892i1,-0.0145 — 0.01401,0.1582 + 0.1195i1,-0.0099 — 

0.00461... 

-0.1943 + 0.0169i1,-0.0331 + 0.07191,-0.0051 —- 0.02531,0.0860 + 
On L501 awa 

0.0135 + 0.0133i1,-0.3556 —- 0.17171,-0.1981 —- 0.0080i1,0.0422 — 

O.OLIAL, a3 

-0.1171 + 0.10941,-0.0243 + 0.03991,0.0955 —- 0.2904i1,0.0060 — 

O12 7 OL y oes 

0.0074 + 0.0105i1,-0.0298 
On LI 88iyees 

0.0298 — 0.1207i1,-0.0030 
0.02251, .<¢ 

-0.0780 — 0.1439i1,-0.1135 —- 0.16401,0.0244 + 0.02831,-0.0426 — 
0.0429i,... 

-0.0690 — 0.06901,0.1682 + 0.19921,0.1395 + 0.23351,-0.0107 
029335 wae 

-0069 + 0.03251,0.0148 + 0.0936i1,-0.0417 —- 0.1614i,-0. 
MOO OSs es&. i 

-0096 — 0.00031,-0.0060 
26405 pees 

-3047 -— 0.23071,-0.0023 
OO31T2; 2% 4 

-0.0375 — 0.01041,0.0088 
O.<O2061; ws % 





+ 0.0183i,-0.0173 + 0.02911,0.0652 — 





+ 0.0140i1,0.0056 —- 0.0015i1,0.0111 + 


0919 


+ 0.0033i1,-0.0985 + 0.07481,0.2024 





+ 0.0014i,-0.0111 + 0.00441,0.0511 
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+ 0.0020i1,0.0043 + 0.0039i,-0.0177 


80 


0.0699 + 0.08841,-0.1185 - 0.1416i1,-0.3045 - 0.29721,-0.0372 - 
0.023217. 
0 
0 





-0511 + 0.00641,-0.0678 + 0.04471,0.0417 - 0.1090i1,-0.0145 —- 





A099, ak 
-0.0286 - 0.07671,0.0594 + 0.06941,-0.0872 - 0.03681,0.0862 - 


0.01731, a4 
0.2500 - 0.2500i1,0.0322 - 0.09931,0.0091 + 0.10871,-0.0406 - 
OOF 9ST piaa a 

0.0542 + 0.04391,0.3314 + 0.07931,0.1083 - 0.02251,-0.0699 - 
0 

0 
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C4921 ems 

.0226 - 0.03521,-0.0035 + 0.02091,0.0373 + 0.11161,0.0257 4 

~Q130L 4.64. 
-0.0104 + 0.01361,-0.0033 - 0.00941,0.0052 + 0.00401,0.1601 + 
O80 3:83 iy vaca ob 
0.1071 - 0.01291,-0.0704 + 0.03171,0.0095 - 0.00751,0.0014 — 
0. O0LS2, «64 
-0.1627 + 0.19901,-0.1191 + 0.12931,0.0786 - 0.06421,0.0073 - 
0.00431, a8 
-0.0194 + 0.00871,0.3171 - 0.12231,0.2274 - 0.09101,-0.0751 + 
0.0409i,... 
—0.0111 + 0.01411,-0.0036 — 0.01501,0.1417 + 0.15831,0.2027 + 
Ou TB 3215 6 « 
-0.0819 - 0.08191,-0.0446 —- 0.06051,0.0189 + 0.04451,-0.0329 























0:26 301 7 aé & 

0.0508 - 0.28561,-0.0300 + 0.0631i1,-0.0340 + 0.03851,0.0112 —- 

0: 00301, a4 

0.0015 - 0.06611,0.0906 —- 0.16621,-0.0337 + 0.04111,-0.0725 + 

0.065615 a08 

0.0271 - 0.01651,-0.0871 + 0.0116i1,-0.0527 - 0.08691,-0.0046 + 
0%. 02191 5.4.48 

-0.0443 + 0.09191,0.0172 - 0.02711,-0.0717 + 0.09851,-0.0956 + 
0.12661; ws « 

0.0018 - 0.00321,-0.0076 - 0.0156i1,0.0044 + 0.0034i1,-0.0200 - 

003491 ye 4's 

-0.0086 — 0.15101,0.0021 —- 0.00821,-0.0705 + 0.14811,0.0110 - 

0. 01591, 464 


-0.0990 + 0.10291,-0.3336 + 0.23601,-0.0441 + 0.01781,0.2044 
0.02581, au 

0.0000,0.1005 + 0.02851,0.3725 + 0.14231,0.0798 + 0.03361,... 
-0.2165 - 0.09451,-0.0218 - 0.01031,-0.0670 - 0.03811,-0.2929 — 


OL2 231 Ty au 
-0.0748 - 0.08221,0.1079 + 0.18091,0.0141 + 0.03971,0.0057 + 
0.036817 48 
-0.0129 + 0.18431,-0.0204 + 0.03031,0.0395 + 0.03821,0.0109 + 
0.0408i,... 
0.0013 + 0.02481,-0.0377 + 0.3309i1,-0.0533 + 0.17671,0.1048 — 
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0.0641 - 0.05781,0.0092 —- 0.00311,0.1928 + 0.02351,0.0948 + 0.0611i,... 


-0.0407 —- 0.08061,0.0165 - 0.04611,0.0010 —- 0.00051,0.2000 + 
Oe OO 222, aa 3 

0.1978 + 0.04821,-0.1638 - 0.0659i,-0.1127 - 0.05971,0.0018 + 
0. 00121 sca 

-0.0248 - 0.02481,0.0599 + 0.03581,-0.0880 - 0.06431,-0.1278 - 
0:0 LO OS ayaa 
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.0041 + 0.00331,-0.1568 - 0.1162i1,-0.1636 - 0.09071,0.0742 + 


.0886 - 0.03771,-0.0008 + 0.00081,0.0790 —- 0.08931,0.1306 - 


-0.0626 + 0.02641,-0.1228 - 0.01871,-0.0029 - 0.0024i1,-0.0642 - 


OES Ty a4 


-0.0783 - 0.28221,0.0120 + 0.10221,-0.0010 + 0.27251,-0.0020 + 


0: OL TS aw a 


-0.0257 + 0.10271,-0.0794 + 0.17001,0.0223 - 0.02481,0.1218 - 


0.0640i,... 





0.0198 - 0.00401,0.0720 —- 0.00041,0.2306 + 0.03771,-0.0258 — 


Os O09 TLS wid 6 


-0.2046 - 0.14531,-0.0303 - 0.04121,-0.0242 - 0.0758i,-0.0144 — 


0:32:60 OL pws 


0 + 0.0000i1,-0.0803 + 0.3119i1,-0.0190 + 0.06731,-0.0037 + 0.03181,... 
0.0604 + 0.08291,0.0173 + 0.00231,-0.2430 + 0.0084i1,-0.1050 — 


0:00 LAS ee % 


-0.0369 —- 0.00621,-0.2672 —- 0.10951,-0.0531 - 0.03941,0.1920 + 


O21 0315. «a0 


0.0670 + 0.09151,0.0081 + 0.00891,0.1187 + 0.03571,0.0513 - 0.0255i,... 


-0.1787 + 0.19121,-0.0991 + 0.14991,-0.0029 + 0.0060i,-0.0600 + 











On. 21461; sie % 

0.0091 + 0.07841,-0.1196 - 0.0696i1,-0.1351 + 0.03411,0.0018 - 
Oi O01 2: pare 

-0.1648 + 0.17441,-0.0797 + 0.14291,0.0440 - 0.26061,-0.0547 - 
0.19490, a4 

0.0032 + 0.00401,-0.1447 - 0.1054i1,-0.1309 - 0.07161,0.1392 + 
OOF ADT yaa a 

0.0507 + 0.05071,0.0008 —- 0.00051,-0.0924 + 0.0075i1,-0.1597 — 
0.00551 5.4-4,6 

0.2167 + 0.02471,0.2508 + 0.05231,0.0017 + 0.00061,0.0670 + 0.0444i,... 
0.0588 + 0.08581,-0.0197 - 0.1036i1,0.0143 - 0.18581,0.0021 —- 
0.00961 sy <4 

0.0272 - 0.08201,0.0851 —- 0.18111,-0.0987 + 0.1418i1,-0.2467 + 
C2218, wes 


-0.0274 + 0.01321,-0.0791 


0:05 825 we % 
0.2861 + 0.18571,0.0276 


+ 0.02511,0.0131 + 0.01 








0.0077 + 0.04201,0.0180 
OO O22 51 5<-%. 6 
-—0.1333 + 0.19661,0.0670 
0.021615 sie 


—- 0.05491,0.3598 - 0.1 


+ 0.01171,-0.2425 - 0.03281,0.1475 + 


L513 70.0110 0400011). 4 is 
+ 0.2580i1,-0.0074 + 0.0690i1,-0.0077 + 


16041,0.0952 —- 


-0.0000,0.0997 -— 0.03051,-0.0055 + 0.01401,0.0949 + 0.1927i,... 


0.0735 + 0.08921,-0.0145 
0.0046i,... 





- 0.01401,0.1582 + 0.1 


1195i,-0.0099 — 


-0.1943 + 0.01691,-0.0331 + 0.07191,-0.0051 - 0.02531,0.0860 + 


On L50T Ly eas 


0.0135 + 0.01331,-0.3556 - 0.1717i1,-0.1981 - 0.00801,0.0422 — 


On OL TAT 7 ace a 


—0.1171 + 0.10941,-0.0243 + 0.03991,0.0955 - 0.29041,0.0060 - 


Pom Dae ilar a 
.0074 + 0.01051,-0.0298 
pel isnsyueear i 
.0298 - 0.12071,-0.0030 
eQ2A D1 yp eas 
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+ 0.0183i,-0.0173 + 0.02911,0.0652 — 


+ 0.0140i1,0.0056 - 0.0015i1,0.0111 + 


-0.0780 - 0.14391,-0.1135 —- 0.16401,0.0244 + 0.02831,-0.0426 - 











0042975 wns ie 

-0.0690 —- 0.06901,0.1682 + 0.19921,0.1395 + 0.23351,-0.0107 - 
0.029317 edo 

0.0069 + 0.0325i1,0.0148 + 0.09361,-0.0417 - 0.16141,-0.0919 - 
O. 1L0997; a4 

0.0096 - 0.00031,-0.0060 + 0.0033i1,-0.0985 + 0.07481,0.2024 — 
0.16401 5-a08 

0.3047 - 0.23071,-0.0023 + 0.0014i1,-0.0111 + 0.00441,0.0511 - 
O'. 005115 «#6 

-0.0375 - 0.01041,0.0088 + 0.00201,0.0043 + 0.0039i1,-0.0177 - 
0.020615 wis 

0.0699 + 0.08841,-0.1185 - 0.1416i1,-0.3045 - 0.29721,-0.0372 - 
OO 232 iy aide ob 

0.0511 + 0.00641,-0.0678 + 0.04471,0.0417 - 0.1090i1,-0.0145 —- 
0.. 405915. 404 

-0.0286 —- 0.07671,0.0594 + 0.06941,-0.0872 - 0.03681,0.0862 - 
0.0173i]'; 
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64; 
.0924; 
64; 
.00035; 
64; 
.0924; 
04 
.0003 
04 
.0925 
04 
.00035; 
64; 
.0924; 
64; 
.00035; 
64; 
.0924; 
64; 
.00035; 
64; 
.0924; 
64; 
.00035; 
64; 
.0924; 
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0005 


;0.143-0.0133 
;0.002-0.1323 
;-0.013+0.1 
;-0.132+0. 
;0.143-0.01 








;-0.132/ 
;0.143-0.01 
;0.002-0.1323 
50.013+0.1434 
;-0.132+0.002 








7=0.1324 
70.143-0.01 








.132+0.0024;-0.013-0.0794;0. 
;-0.013-0.0794;-0. 
Mchkc sre 
3002-051 
.143-0. 
.132+0. 
.013+0.1 
.002-0.1 
.143-0. 
AB: 
.013+0.1 
.002-0.1 
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Se Ge ee 
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se Ue ute 





-0.079-0.01335 
7-0.079-0.013 
;7-0.013-0.079 
-0.013-0.0795 
-0.079-0.0133 
7-0.079-0.013 
;-0.013-0.079 
-0.013-0.0795 
-0.079-0.01335 
7-0.079-0.013 
;7-0.013-0.079 
-0.013-0.0795 
-0.079-0.01335 
-0.079-0.0133 
7-0.013-0.079 
-0.013-0.0795 
-0.079-0.01335 
7-0.079-0.013 
7-0.013-0.079 
-0.013-0.0795 
-0.079-0.01335 
7;-0.079-0.013 
7-0.013-0.079 
-0.013-0.07935 
-0.079-0.0133 
7-0.079-0.013 
7-0.013-0.079 
-0.013-0.0795 
-0.079-0.01335 
7-0.079-0.013 
7-0.013-0.079 
-0.013-0.07935 
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0.046+0.0464;0.002-0.13245;-0.079-0.0135;-0.013+0.143); 
0.000+0.0924;-0.013+0.1435;-0.079-0.0135;0.002+0.132); 
0.046+0.0464;-0.132+0.00235;-0.013+0.0795;0.143-0.013); 
0.092+0.0005;0.143-0.01345;-0.013-0.0795;-0.132+0.0023; 
0.046+0.0464;0.002-0.13245;-0.079-0.0135;-0.013+0.143); 
0.000+0.0924;-0.013+0.1435;-0.079-0.0135;0.002-0.132); 
-—0.078+0.0003;0.012-0.0984;0.092-0.1064;-0.092-0.115j; 
-—0.003+0.0545;0.075+0.0744;-0.127+0.02145;-0.122+0.0174; 
—0.035+0.1515;-0.056+0.02275;-0.060-0.08145;0.070-0.014j; 
0.082-0.0924;-0.131-0.0653;-0.057-0.03945;0.037-0.098); 
0.062+0.06235;0.119+0.0045;-0.022-0.1615;0.059+0.0153;. 
0.024+0.0594;-0.137+0.0475;0.001+0.1155;0.053-0.0043;. 
0.098+0.0264;-0.038+0.1063;-0.115+0.0555;0.060+0.088); 
0.021-0.0283;0.097-0.0835;0.040+0.1115;-0.005+0.1203;. 
0.156+0.0003;-0.005-0.1205;0.040-0.11135;0.097+0.0833;. 
0.021+0.0284;0.060-0.0884;-0.115-0.0554;-0.038-0.106); 
0.098-0.0263;0.053+0.0045;0.001-0.11535;-0.137+0.0473;. 
0.024-0.0593;0.059-0.0154;-0.022+0.1615;0.119-0.0047;... 
0.062-0.06235;0.037+0.09845;-0.057+0.0395;-0.131+0.065;... 
0.082+0.09245;0.070+0.01445;-0.060+0.0815;-0.056-0.022)4; 
-—0.035-0.1515;-0.122-0.0174;-0.127-0.02145;0.075-0.0743; 
-—0.003+0.0545;-0.092+0.1154;0.092+0.10645;0.012+0.0984; 
-—0.156+0.00053;0.012-0.0984;0.092-0.1064;-0.092-0.115j; 
-—0.003-0.0545;0.075+0.0744;-0.127+0.02145;-0.122+0.01735; 
-—0.035+0.1515;-0.056+0.0224;-0.060-0.08145;0.070-0.0149; 
0.082-0.0924;-0.131-0.0654;-0.057-0.03945;0.037-0.098); 
0.062+0.06235;0.119+0.0045;-0.022-0.16135;0.059+0.0153;. 
0.024+0.05935;-0.137+0.0475;0.001+0.1155;0.053-0.0043;. 
0.098+0.0264;-0.038+0.1063;-0.115+0.05545;0.060+0.088); 
0.021-0.02835;0.097-0.0835;0.040+0.1115;-0.005+0.1203;. 
0.156+0.0003;-0.005-0.1205;0.040-0.11145;0.097+0.0833;. 
0.021+0.0284;0.060-0.0884;-0.115-0.0554;-0.038-0.106); 
0.098-0.0263;0.053+0.0045;0.001-0.11535;-0.137-0.0473;. 
0.024-0.0593;0.059-0.0155;-0.022+0.1615;0.119-0.0043;. 
0.062-0.06245;0.037+0.09845;-0.057+0.03943;-0.131+0.065); 
0.082+0.0924;0.070+0.01435;-0.060+0.0815;-0.056-0.022)3; 
-0.035-0.1515;-0.122-0.0174;-0.127-0.02145;0.075-0.0743; 
—0.003+0.0545;-0.092+0.11545;0.092+0.10645;0.012+0.0984; 
—0.156+0.0005;0.012-0.0985;0.092-0.1063;-0.092-0.115j;... 
-0.003-0.0545;0.075+0.0745;-0.127+0.02135;-0.122+0.01745;... 
—0.035+0.1515;-0.056+0.0225;-0.060-0.0813;0.070-0.014j;... 
0.082-0.0924;-0.131-0.0654;-0.057-0.0395;0.037-0.09874;... 
0.062+0.06235;0.119+0.0045;-0.022-0.1615;0.059+0.0157;... 
0.024+0.0594;-0.137+0.0475;0.001+0.1155;0.053-0.00474;... 
0.098+0.0264;-0.038+0.1065;-0.115+0.0555;0.060+0.08874;... 
0.021-0.0283;0.097-0.0835;0.040+0.1115;-0.005+0.12035;]; 

% WiFi Auto FAM (After [11]) % 
FUNCE LON 

[Sxp, fo,Np, centrow, centcol]=WiFi_autofam(famin,fs,df,dalpha,ss,... 


index,v, shift); 
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if nargin~=8 
error('Wrong number of arguments.'); 
end 
for ww=6:index*2 
xin(1, (ss-16) * (ww-1) +1: (ss-16) *ww) =famin(1,ss* (ww-1)+16+1:ss*ww) ; 
end 
sdefine parameters 
Np=pow2 (nextpow2 (fs/df)); 
L=Np/4; 
P=pow2 (nextpow2 (fs/dalpha/L) ) ; 
N=P*L; 
Sxp=zeros (Np+1,2*N+1); 
for w=1:l:index/2 
x=xin(1, (shift) * (w-1)+1: (shift) *w); 


Sinput channalization 
if length (x) <N 


x (N) =0; 
elseif length(x)>N 
x=x(1:N); 
end 
NN=(P-1) *L+Np; 
XX=X} 
xx (NN) =0; 
XX=XX(:); 


X=zeros (Np,P); 
for k=0:P-1 

X(:,k+1)=xx (k*L+1:k*L+Np) ; 
end 


Swindowing 


a=hamming (Np) ; 
XW=diag (a) *X; 
SXW=X; 


Sfirst FFT 


XF1l=fft (XW); 
XFl=fftshift (XF1); 
XF=[XF1(:,P/2+1:P) XF1(:,1:P/2)]; 


Sdownconvertions 
E=zeros (Np,P); 
for k=-Np/2:Np/2-1 
for k0=0:P-1 
EB (k+Np/2+1,k0+1)=exp (-i*2*pi*k*kO*L/Np) ; 











end 
XD=XF1.*E; 
XD=conj (XD'); 





smultiplication 
XM=zeros(P,Np%2); 
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for k=1:Np 
for c=1:Np 
XM(:, (kK-1) *Npt+c)=(XD(:,k) .*conj(XD(:,c))); 
end 
end 


Ssecond FFT 
XF2=fft (XM); 
XF2=fftshift (XF2); 
XF2=[XF2(:,Np*2/2+1:Np*2) XF2(:,1:Np*2/2)]; 
XF2=XF2 (P/4:3*P/4,:); 
M=abs (XF2); 
alphao=-fs:fs/N:fs; 
fo=-fs/2:fs/Np:fs/2; 
Sx=zeros (Np+1,2*N+1); 
for k1l=1:P/2+1 
for k2=1:Np%2 
if rem(k2,Np) == 
c=Np/2-1; 
else 
c=rem(k2,Np) -Np/2-1; 
end 
k=ceil (k2/Np) -Np/2-1; 
p=k1-P/4-1; 
alpha=(k-c) /Npt+ (p-1) /L/P; 
f=(kt+c) /2/Np; 
if alpha<-1 | alpha>l 
k2=k2+1; 
elseif f<-.5 | f>.5 
k2=k2+1; 

















else 
kk=1+Np* (f+.5); 
11=1+N* (alphat+l) ; 
Sx (round (kk), round(11))=M(k1,k2); 
end 
end 
end 
SXp=SxXptSx; 
end 
Sxp=Sxp./max (max (Sxp) ); 





oe 


Surface Plot 


Alpha=num2str (round (1+N* (alpha/fs+1))); 
Fo=num2str (floor ((Npt+1)/2)+1); 
V=num2str (v); 

figure (4) 

surfl(alphao, fo,Sxp); 

view (-37.5, 60); 

















xlabel('alpha'); 
ylabel('f'); 
zlabel('Sx'); 


colormap jet 
colorbar 
axis([0 20*10*%6 —9*10%6 9*10%6 O 1]) 
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SContour Plot 


figure (5) 
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bar 





colormap winter 


color 


SCross-Section Plot 


figure (6) 


)); 


t (alphao, Sxp (floor ( (Npt+1) /2) +1 


, 


round (1+N* (alpha/fs+1)))) 
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Baseband Data Setup and Demodulator 
Written by Steven Schnur 


7: 


(L,N,m,Ndata,CP,y,bk, Frame) 


WiFi_BasebandDemod 


[ck] 


Function 


if nargin~=8 


error('Wrong number of arguments") 


end 
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Remove Preamble 
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:N+1],Frame)... 


& mod([-2 


:N],Frame) 


=mod ( [-3 


index 


sCreats logic 4preamble 


N+3],Frame) ; 


& mod([0 


:N+2],Frame) 


& mod([-l 


insertion 


pad 


Loop to break data string into L-fft sized 


% 


for b 


blocks 


Loop skips demod code if preamble logic 


% 


satisfied 


ole 


else if 


pad=padt+1; 


Loads L# I-Q data into vector 


% 
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Removes CP 
L-point FFT converts time samples to I-Q data 


6 


fe) 


% 


yl (CP+1:L+CP,1); 


YR=fft (y2,L); 


y2 


oe 


Extracts user data from appropriate SC 


ol? 


l, 





ol CY ame 
ef 


"hard decision'); 





aLCy aem 





“CG 
(3+33) 
(=l=19) 


(1+13) 








noise 
(=3=35)) 


of Cry 
sof Cry 


Plot recieved I-Q data with 
fC; 


NdataxNxm bits column wise 


[ 


(-3+3)) 


"hard decision'); 


et ae eS  . 


(3-13) 





Normalization Value of 16QAM Sym 





° 
Oo 


of Cy 
if: Cy (H14195) 


sf Cy 
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(1-33) 
(-3-1)) 


(3417) 





DecisionType', 
"DecisionType', 


~~~ He 


af Cy 
if Cy 
eae; 


Generates 640AM Demodulation Object 
fC 








Generates 16QAM Demodulation Object 


(1+33) 
(~3+19) 





(-1-3)) 





41.65.51) 


yer rewrvr wv wr bt 


ef Cy 
vf Cy 


YR(23227,1)7 
"OutputType', 


YR(9221,.1) 7 


af-Cy 


ee 


modem.gengamdemod('Constellation', 


)=YR(45:57,1); 


19:24,1)=YR(59:64,1); 


figure (2) 





)=YR(39:43,1); 





("Recieved Quadrature Data") 


=14.5 125 


('Recieved In-Phase Data") 
grid on 


(3-39) 
(-1+34) 
(1-13) 





hrec 


~e rye rer vr ve vee YS YS YS WS 


25:30,1)=YR(2:7,1); 


31:43,1) 





44:48,1) 


~erYye rye ve wv YH 








sqrt (10); 


sqrt (42); 
hrec=modem.gengamdemod('Constellation', 





title(['Recieved I-Q Data']) 


plot. (y¥1, **"); 
xlabel 

ylabel 

axis ( 

if m==6 

else if m==4 
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else if m==2 


sqrt (2); 


(1-14) /c, 


(-1-15) /¢; 


[ 





hrec=modem.genqamdemod('Constellation', 


"DecisionType', ‘hard decision'); 


"Bat, 


"OutputType', 
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Demodulates received data to baseband IQ data 
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End of for bl 
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end 


m*pad*Ndata) =Yp; 


ck (1,m*Ndata* (pad-1) +1 


=Y1p; 


Y2 ((pad-1) *Ndata+1l:pad*Ndata,1) 


end 
end 


0) loop 


(index (1,b) 


End of if 


% 


N loop 


=1: 


End: of tor 6 


% 


a 


length (ck) ) 


bk=bk (1,1 


Compares rec data string to trans 
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max (cumsum(err) ) / (m*Ndata*N) 


num2str(B 
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figure (3) 
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NdataxNxm bits column wise 
Plot recieved I-Q data with noise 
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Baseband Data Setup and Modulator 
Written by Steven Schnur 


oe 
AP Al? ol? 


(L,N,m, Ndata,bk,CP,G, Frame) 





WiMax _BasebandMod 


[x] 


Function 


if nargin~=8 


error('Wrong number of arguments") 


end 


Generates 64Q0AM Modulation Object 


if m==6 


c 





Normalization Value of 64QAM Sym 


% 


(3+35)/c, (3+4)/c,... 


[ 


modem.gengammod('Constellation', 


object 
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Generates 16QAM Modulation Object 
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else if m==4 


Normalization Value of 16QAM Sym 
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sqrt (10) 
object = modem.gengqammod('Constellation', 
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else if m==2 


sqrt (2) 


object 


(-1+3) /G; 
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IQ data 


Ndata loop 


ub 


to load m-by-Ndata fft vector 


Grabs one m-symbol worth of data 
ding of storage vect with 1 m-sym 


Modulates m sized sym with obj 


Loop 
Loa 
Stacks IQ data in a column vector 
End of w 
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% 
% 
% 
% 


% 


, 


Ndata 


il. 
modulate (object,mk') 


Xp 


for w 
end 
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% Pilot Subcarrier Sequence Generator % 
PilotSeqGen =seqgen.pn('GenPoly', 

[10000000010 1],... 

‘Iiratialstates’, [1 111 2) 112 1 2 1 1) pees 

"NumBitsOut',1); 
set (PilotSeqGen, 'NumBitsOut', N); 
wk=generate (PilotSeqGen) ; % Creates PN Seq for Pilot SubCar 
for k1=1:N; & Symbols as WIMAX Standard 
if wk(k1,1)==0; % Logic loop to create compliment 

wkcomp (k1,1)=1; % of PN Seq for Negitive Pilot SC 
else if wk(k1,1)==1; 

wkcomp (k1,1)=0; 

end 
end 
Negpilot (k1,1)=1-2*wk (k1,1); % CreatesWIMAX Standard Pilot SC 
Pospilot (k1,1)=1-2*wkcomp (k1,1); % for pos and neg Pilot SC 
end 
NPSCsym=Negpilot (k,1)j; % Loads Pilot data sym by sym 
PPSCsym=Pospilot (k,1); 
% Load Vector X with Guard SC=0, DC=0, Pilot SC and Data in % 
% appropriate SC locations % 









































X(1,1)=0; % DC Null 
X(2:12,1)=X1(97:107,1); % Data 
X(13,1)=G*PPSCsym; % Positive Pilot SC 
X (14:37, 1)=X1(108:131,1); % Data 
X(38,1)=G*PPSCsym; & Positive Pilot SC 
X(39:62,1)=X1(132:155,1); % Data 
X(63,1)=G*PPSCsym; % Positive Pilot SC 
X (64:87,1)=X1(156:179,1); % Data 
X(88,1)=G*PPSCsym; % Positive Pilot SC 
X(89:101,1)=X1(180:192,1); % Data 
X(102:128,1)=0; % Lower Guard Band 
X(129:156,1)=0; % Upper Guard Band 
X(157:168,1)=X1(1:12,1); % Data 
X(169,1)=G*NPSCsym; % Negitive Pilot SC 
X(170:193,1)=X1(13:36,1); % Data 
X(194,1)=G*NPSCsym; %& Negitive Pilot SC 
X(195:218,1)=X1(37:60,1); % Data 
X(219,1)=G*NPSCsym; %& Negitive Pilot SC 
X(220:243,1)=X1(61:84,1); % Data 
X(244,1)=G*NPSCsym; % Negitive Pilot SC 
X (245:256,1) =X1(85:96,1); & Data 
Xplot ((k-1) *L+1:k*L,1) =X; % Dummy Vector for plotting 
end & End of else if index(1,k)==1 loop 
end 6 End if ((index(1,k)==0) & (index(1,k+1)==1)) loop 
end 6 End if ((index(1,k)==0) & (index(1,k+1)==0)) loop 
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io) 
% IFFT Baseband Modulated data to generate time dat 
SESESESE SEES EES SE EEE EES EEEEE SEES EEEE ESE EEEESE ESC EEEEESEEEEEEEESEEEEESSESEESEESS 








xt=ifft (X,L); & L size IFFT to convert I-Q data to time samples 
g(1:CP,1)=xt ((L-CP+1):L,1); % Loads CP worth of Data from end of sym 
x1(1:CP,1)=9; % Loads CP to front of vector 

x1 (CP+1:L+CP,1)=xt; % Loads rest of OFDM Symbol into vector 

x (((k-1) * (L+CP) +1) :k* (L+CP),1)=x1l; % Fills transmit vector w/CP app 
end % End for k=1:N loop 


figure (1) 

plot (xX, 7%"); 6 Plot Tx Constellation%Stitle('Transmitted I-Q Data") 
xlabel('Transmitter In Phase Data') 

ylabel('Transmitted Quadrature Data') 

axis([=125 225. =1.5 155] )¢ 

grid on 














CSCSSCCSEEEEESEEEEEEEEEEEEEEEE SESE SESS SSS SSS SSC SESEEEEEEEEEEEEEEEESESESSESSES 
ie) coy 
% WiMax Preamble Construction: First Symbol % 
o9O0 aoo00000000 oo0 ago00g0000cgcCoCoOoCoOoOoOoOOOOoOOOOWOOoOoOoOOO 8 gooog00g000ccogccCcCoOoCoCoCACNNND 


function [X]=preamble_1_64; 


pall=[ 1 - 11,1 - 1i,-1 - 11,1 + 11,1 - 11,1 - 11,-1 + 11,1 - 11,1 - 11,1 - 11,1 + 11,-1 - 

1i,1 + 11,1 + 11,-1 - 11,1 + 11,-1 - 14,-1 - 11,1 - 14,-1 + 11,1 - 14,1 - 11,-1 - 11,1 + 1i,1 - 11,1 - 
1i,-1 + 1i,1 - 1i,1 - 11,1 - 14,1 + 14,-1 - 14,1 + 11,1 + 1i,-1 - 14,1 + 14,-1 - 11,-1 - 11,1 - 14,-1 
+ ijl - 1i,1 - 1i,-1 - 14,1 + 11,1 - 14,1 - 14,-1 + 14,1 - 14,1 - 11,1 - 11,1 + 14-1 - 14,1 + 1i,1+ 
hi-t = 14,1 diel = 14,-1°> Ia,1 = 11,-1 + 11,1 11,1 la] = 1aj-1-+ 1,1 + 14,1 + 11,-1 = 14,1 
+ 1i,] + 11,1. 1ij-1 + 11,1 = Lij=1 + 14,-1 + 11,1 = 14,-1 + 11,1 = 11,1 = 11,1 + 11,-1 = 11,-1 = 
1i,-1 - 14,-1 + 11,1 - 1i,-1 - 14,-1 - 11,1 + 11,-1 - 14,-1 - 11-1 - 11,1 - 14,-1 + 11,1 - 14,1 - 11,-1 
+ 1i,1 - 1i,-1 + 14,-1 + 14,-1 - 11,1 + 11,0 + 01,-1 - 11,1 + 14,-1 + 1i,-1 + 1,-1 - 11,1 + 11,1 + 
1i,1 + 1i,-1 - 11,1 + 14,1 - 14,1 - 14,1 - 11,-1 + 14,-1 + 11,-1 + 11,-1 + 11,1 - 14,-1 - 114,-1 - 11,- 
1+ 1i,1 - 1i,1 + 11,1 + 11,-1 + 11,1 - 11,1 - 11,1 - 14,-1 + 14,1 - 11,-1 - 14,-1 - 1i,-1 - 11,1 + 
1i,1 + 1i,1 + 14,1 + 14,-1 - 14,-1 + 1i,-1 + 14,1 + 11,-1 - 11,1 - 14,1 - 14,1 + 14,-1 - 114,-1 - 11,- 
1 - 14,1 + 1i,-1 - 14,-1 + 11,-1 + 11,-1 + 1i,1 - 11,1 - 14,1 - 11,1 - 1i,-1 + 11,1 + 11,1 + 11,-1 - 
1i,1 + 1i,-1 + 11,-1 + 14,-1 - 11,1 + 11,1 + 11,1 + 11,-1 - 11,1 + 11,1 - 11,1 - 11,1 - 11,-1 + 11,- 
1+ 1i,-1 + 1i,-1 + 11,1 - 14,-1 - 14,-1 - 14,1 - 1i,-1 + 11,-1 - 11,-1 - 11,1 - 14,-1 + 14-1 + 11,-1 
lil = Thal lad + il a a el = Tel = lel = Tiel = Ia Tis ta = Ta; 


n=length (pall); 
a=0; 
ind=mod([0:n],4); 
for k=1:n 
if ind(k) == 
a=atl; 
pinsert51 (a,1)=sqrt (2) *conj (pall (k,1)); 
end 
end 
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p64=[pinsert51(1:25,1)',pinsert51(27:51,1)']; 
Spreamble4_64=[zeros(1,102:128),p64,p64,p64, p64]; 
SX=preamble4_64'; 

X(1,1)=0; 

X(2:51,1)=p64; 

X(52:101,1)=p64; 

X(102:128,1)=0; 

X(129:156,1)=0; 

X(157:206,1) =p64; 

X(207:256,1) =p64; 





2-9-9. 2. 9. 9.2. 9. 9.'9- 


function [X]=preamble_2_128; 

clear 

cle 

pall=[1 - li,1 - li,-1 - 14,1 + 11,1 - 14,1 - 1i,-1 + 14,1 - 11,1 - 11,1 - 11,1 + 11,-1 - 

1i,1 + 1i,1 + 1i,-1 - 1i,1 + 11,-1 - 1i,-1 - 11,1 - 1i,-1 + 14,1 - 1i,1 - 11,-1 - 11,1 + 1i,1 - 11,1 - 
1ij-1 + 1i,1 - 14,1 - 14,1 - 14,1 + 114,-1 - 14,1 + 14,1 + 14,-1 - 11,1 + 114,-1 - 14,-1 - 11,1 - 11,-1 
+ 14,1 - 11,1 - 14,-1 - 14,1 + 1i,1 - 14,1 - 1i,-1 + 11,1 - 14,1 - 14,1 - 14,1 + 14,-1 - 14,1 + 1i,1 + 
1ij-1 - 11,1 + 114,-1 - 1i-1 - 11,1 - 1i,-1 + 14,1 + 14,1 + 14,1 - 11,-1 + 14,1 + 11,1 + 114,-1 - 11,1 
+ 1i,1 + 1i,1 + 1i,-1 + 1i,1 - 11,-1 + 1i,-1 + 11,1 - 1i,-1 + 11,1 - 14,1 - 11,1 + 11,-1 - 11-1 - 
1i,-1 - 1i,-1 + 1i,1 - 1i,-1 - 11,-1 - 14,1 + 11,-1 - 1i,-1 - 1i,-1 - 14,1 - 1i,-1 + 14,1 - 11,1 - 11,-1 
+ 1i,1 - 11,-1 + 11,-1 + 114,-1 - 11,1 + 11,0 + O1,-1 - 11,1 + 1i,-1 + 11,-1 + 14,-1 - 11,1 + 11,1 + 
lil + 1i,-1 - 1ij1 + 11,1 - 1i,1 - 14,1 - 1i,-1 + 1i,-1 + 1i,-1 + 11,-1 + 11,1 - 1i,-1 - 11,-1 - 11,- 
1+ 1i,1 - 11,1 + 11,1 + 1i,-1 + 1a,1 - 11,1 - 11,1 - 14,-1 + 11,1 - 14-1 - 1i,-1 - 11,-1 - 1i,1 + 
1i,) + 11,1 + 11,1 + 1i,-1 - 1i,-1 + 1i,-1 + 14,1 + 1i,-1 - 14,1 - 11,1 - 11,1 + 11,-1 - 1i,-1 - 1,- 
1 - 1i,1 + 1i,-1 - 1i,-1 + 1i,-1 + 1i,-1 + 1a,1 - 11,1 - 11,1 - 14,1 - 1i,-1 + 14,1 + 11,1 + 11,-1 - 
adel asl els i el day ae a a at = Ta = as te ae 
1+ 1i,-1 + 11,-1 + 1i,1 - 1i,-1 - 11,-1 - 1i,1 - 1i,-1 + 1i,-1 - 1i,-1 - 11,1 - 11,-1 + 14,-1 + 11,-1 
+ 1i,d - Lij-1 + 1i,1 + 14,1 + 11,1 + 1i-1 - 1i,-1 - 14-1 - 11,-1 - 14,1 + 11,1 - 14,1 - 11]; 

















n=length (pall); 

a=0; 

ind=mod([0:n],2); 

for k=l1:in 

if ind(k) == 
a=atl; 
pinsert1l01(a,1)=sqrt (2) *conj (pall (k,1)); 

end 

end 


pl28=[pinsert101(1:50,1)', pinsert101(52:101,1)'J; 
spreamble2_128=[p128,p128]; 

sX=preamble2_128'; 

X(1,1)=0; 

X(2:101,1)=p128; 
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X(102:128,1)=0; 
X(129:156,1)=0; 
X(157:256,1)=p128; 


function [Sxp, fo,Np, centrow, centcol]=WiMax_autofamtest (famin,fs,df,... 
dalpha,ss,index,CPl1, shift); 

if nargin~=8 
error ('Wrong number of arguments.') 

end 


for ww=3:index*2 
xin(1, (ss—-CP1) * (ww-1)+1: (ss-CP1) *ww) =famin(1,ss* (ww- 
1)+CP1+1:ss*ww); 
end 
sdefine parameters 
Np=pow2 (nextpow2 (fs/df)); 
L=Np/4; 
P=pow2 (nextpow2 (fs/dalpha/L) ) ; 
N=P*L; 
Sxp=zeros (Np+1,2*N+1); 
for w=1:1:floor(index/2) 
x=xin(1, (shift) * (w-1)+1: (shift) *w); 


Sinput channalization 

if length (x) <N 
x (N) =0; 

elseif length(x)>N 
x=x(1:N); 

end 

NN=(P-1) *L+Np; 

XX=X} 

xx (NN) =0; 

XX=XX(:); 

X=zeros (Np,P); 

for k=0:P-1 
X(:,k+1)=xx (k*L+1:k*L+Np) ; 

end 


Swindowing 


a=hamming (Np) ; 
XW=diag (a) *X; 


Sfirst FFT 


XF1l=fft (XW); 
XFl=fftshift (XF1); 
XF=[XF1(:,P/2+1:P) XF1(:,1:P/2)]; 
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Sdownconvertions 
E=zeros (Np,P); 
for k=-Np/2:Np/2-1 
for k0=0:P-1 
EB (k+Np/2+1,k0+1)=exp (-i*2*pi*k*kO*L/Np) ; 











end 
end 
XD=XF1.*E; 
XD=conj (XD'); 





Ssmultiplication 

XM=zeros (P,Np*2); 

for k=1:Np 
for c=1:Np 

XM(:, (k-1) *Npt+c)=(XD(:,k).*conj(XD(:,c))); 





end 
end 


Ssecond FFT 
XF2=fft (XM); 
XF2=fftshift (XF2); 
XF2=[XF2(:,Np*2/2+1:Np*2) XF2(:,1:Np%2/2) ]; 
XF2=XF2 (P/4:3*P/4,:); 
M=abs (XF2); 
alphao=-fs:fs/N:fs; 
fo=-fs/2:fs/Np:fs/2; 
Sx=zeros (Np+1,2*N+1); 
for k1=1:P/2+1 
for k2=1:Np%2 
if rem(k2,Np) == 
c=Np/2-1; 
else 
c=rem(k2,Np) -Np/2-1; 
end 
k=ceil (k2/Np) -Np/2-1; 
p=k1-P/4-1; 
alpha=(k-c) /Npt+ (p-1) /L/P; 
f=(kt+c) /2/Np; 
if alpha<-l1 | alpha>l 
k2=k2+1; 
elseif f<-.5 | f>.5 
k2=k2+1; 

















else 

kk=1+Np* (f+.5); 

11=1+N* (alphat+l); 

Sx (round (kk), round(11))=M(k1,k2); 
end 





end 
end 
SXp=SxXptSx; 
end 


Sxp=Sxp./max (max (Sxp) ); 
centrow=floor((Npt1)/2) +2; 
centcol=round(1+N* (alpha/fst1)); 
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2S. 'S- 2 2. ¢ 9. 9. 9. 0-0-9. 2-2. 9.9.9. 2-9. 


FuNnceLon 
[Sxp, fo,Np, centrow, centcol]=WiMax_autofam(famin,fs,df,dalpha,... 

ss,index,v,shift,CP2); 
if nargin~=9 

error ('Wrong number of arguments.') 
end 
for ww=3:index*2 

xin(1, (ss-CP2) * (ww-1)+1: (ss-CP2) *ww) =famin(1,ss* (ww- 
1)+CP2+1:ss*ww) ; 
end 
sdefine parameters 
Np=pow2 (nextpow2 (fs/df)); 
L=Np/4; 
P=pow2 (nextpow2 (fs/dalpha/L) ) ; 
N=P*L; 
Sxp=zeros (Np+1,2*N+1); 
for w=1:1:floor (index/2) 

x=xin(1, (shift) * (w-1)+1: (shift) *w); 


Sinput channalization 
if length (x) <N 
x (N) =0; 
elseif length(x)>N 
x=x(1:N); 
end 
NN=(P-1) *L+Np; 


X=zeros (Np,P); 
for k=0:P-1 

X(:,k+1)=xx (k*L+1:k*L+Np) ; 
end 


Swindowing 


a=hamming (Np) ; 
XW=diag (a) *X; 


Sfirst FFT 


XF1=ff£t (XW) ; 
XFl=fftshift (XF1); 
XF=[XF1(:,P/2+1:P) XF1(:,1:P/2)]; 


Sdownconvertions 
E=zeros (Np,P); 
for k=-Np/2:Np/2-1 
for k0=0:P-1 
EB (k+Np/2+1,k0+1)=exp (-i*2*pi*k*kO*L/Np) ; 
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end 
end 
XD=XF1.*E; 
XD=conj (XD'); 





smultiplication 

XM=zeros (P,Np*2); 

for k=1:Np 
for c=1:Np 

XM(:, (k-1) *Npt+c)=(XD(:,k).*conj(XD(:,c))); 





end 
end 


Ssecond FFT 
XF2=fft (XM); 
XF2=fftshift (XF2) ; 
XF2=[XF2(:,Np*2/2+1:Np*2) XF2(:,1:Np%2/2) ]; 
XF2=XF2 (P/4:3*P/4,:); 
M=abs (XF2); 
alphao=-fs:fs/N:fs; 
fo=-fs/2:fs/Np:fs/2; 
Sx=zeros (Np+1,2*N+1); 
for k1=1:P/2+1 
for k2=1:Np%2 
if rem(k2,Np) == 
c=Np/2-1; 
else 
c=rem(k2,Np) -Np/2-1; 
end 
k=ceil (k2/Np) -Np/2-1; 
p=k1-P/4-1; 
alpha=(k-c) /Npt+(p-1) /L/P; 
f=(kt+c) /2/Np; 
if alpha<-l1 | alpha>l 
k2=k2+1; 
elseif f<-.5 | f>.5 
k2=k2+1; 

















else 

kk=1+Np* (f+.5); 

11=1+N* (alphat+l); 

Sx (round (kk), round(11))=M(k1,k2); 
end 





end 
end 
SxXp=SxXptSx; 
end 


Sxp=Sxp./max (max (Sxp) ); 
% Surface Plot 
centrow=floor((Npt+1)/2) +2; 
centcol=round(1+N* (alpha/fst1)); 
Alpha=num2str(centcol); 
Fo=num2str(centrow) ; 
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figure (4) 


¥ 


surfl(alphao, fo, Sxp) 


view (-37.5, 60) 


¥ 


i 


=2*10%6 2*10%6 0 17) 


0-3.5* 106 





bel ('alpha'); 














bar 





colormap je 
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axis ( 


sContour Plot 
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color 
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(L,N,m,Ndata,CP,y,bk, Frame) 


WiMax_BasebandDemod 


[ck]= 


FUNnCELOTN 


if nargin~=8 


error('Wrong number of arguments") 


end 
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Logic vector to determine when to remove 
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:N],Frame) 


=mod ([-1 


index 


Preamble 


le 


pad 


Loop to break data string into L-fft sized 


% 


for b 


blocks 


Loop skips demod code if preamble logic 


% 


satisfied 


else if 


pad=padt1; 
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yl=y ((b-1) * (L+CP) +1:b* (L+CP),1); % Loads L# I-Q data into vector 
y2=yl(CP+1:L+CP,1); % Removes CP 


— . fo) _ : ‘ - 
YR=fft (y2,L); % L-point FFT converts time samples to I-Q data 
COSC CCC SE EEEEEEEEEEEEEEEEESEE EEE SESS SSS SS SSS SESE SESE EEESESEEEEEEEEEEESES 
co} °o 
% Extracts user data from appropriate SC % 
CESS SCCCSEEEEEEEEEEEEEEEEEEEEE SESE SESS SSS SSS SSS SESEEEEEEEEEEEEEEEEESESESSES 


97:107,1)=YR(2:12,1); 
108:131,1)=YR(14:37,1) 
1323155,1)=YR(39262,-1) 
156:179,1)=YR(64:87,1) 
180:192,1)=YR(89:101,1 
1:12,1)=YR(157:168,1); 
13:36,1)=YR(170:193,1) 
37:60,1)=YR(195:218,1); 
61:84,1)=YR(220:243,1) 
85:96,1)=YR(245:256,1) 
igure (2) 

lot (Y1,'*"); 
itle(['Recieved I-Q Data']) 

label ('Recieved In-Phase Data") 
abel ('Recieved Quadrature Data') 
axis: ([—_L.5 1.5 =1.5 1.5] ) 

grid on 





, 

NdataxNxm bits column wise 

Plot recieved I-Q data with 
noise 


ol? 
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c=sqrt (42); 
hrec=modem.gengamdemod('Constellation', [ (3+33)/c, (3+j)/c,... 
(S455) /ep (37g) fd, (S=39)/cp (3=j) fe, (3-55) /c, (3-75) /¢;, 
































(1+353)/c, (1+3)/c, (1+53)/c, (1+735)/c, (1-33) /c, (1-3) /c, (1-59) /c, 
(1-75) /c, (5+345)/c, (5+35)/c, (5+545)/c, (5+745)/c, (5-345) /c, (5-5) /c, 
(5-53) /c, (5-753)/e (7+33)/c, (7+35)/c, (7453) /c, (74795) /c, (7-395) /c, 
(Jj) /S,. (759) 4Cr (Cl=79) fe, (= 3437) fey (343) /cp (31576 aus 
(oT) (Gr. (ses) i ey We s=i) ee, (eors i) ce. a1) fe, eo) fe, 
(=Ltj)/e, (=it59)/e, (-l+Ti) fe; (= 139) /er. (41-9) /e, (F595) (Cy axe 
(=1=75) 7G, (S435) /ey (5453) /e;- 5457) /e,. S77) ee, (-5=35) LG 
(=5=9)7¢e, (=o-59)/¢, (-S-7 3) /e; (-7435)./¢); (749) 76, (-1+55) fc; 
(-74745)/ce, (-7-34)/e, (-7-3)/ce, (-7-53)/c, (-7-74)/e 1, 
"OutputType', 'Bit', 'DecisionType', ‘hard decision'); 

% Generates 16Q0AM Demodulation Object % 


2.9. 2. 2. 9.2. 9. 2. 9. 9. 9. 8. S-'2- 9. 8. 9-9. O- G- 


else if m== 
hrec=modem.gengamdemod('Constellation', [ (1+ j)./sqrt(10), 
(1+33)./sqrt (10), (1-j)./sqrt(10), (1-34) ./sqrt(10), 
(3+3) ./sqrt (10), (3+33)./sqrt (10), (3-4) ./sqrt (10), (3-33) ./sqrt (10), 
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(-1-j)./sqrt (10), 
($3439) «/SQLE (10) pews 


(-1+33)./sqrt (10), 


(-1+4)./sqrt (10), 
(-1-33)./sqrt (10), 


(-3+ 3) ./sqrt (10), 


"DecisionType', ‘hard decision'); 


"Bae; 


"OutputType', 


oe 


Generates QPSK Demodulation Object 


le 


else if m==2 





hrec=modem.genqamdemod('Constellation', 
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"hard decision') 


"DecisionType', 





'Bit', 


"OutputType', 


oe 


Generates BPSK Demodulation Object 
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else if m==1 
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"hard decision') 





hrec=modem.genqamdemod('Constellation', 
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"DecisionType', 


"Blt, 


‘OutputT ype", 
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End of else if 


% 
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Demodulates received data to baseband IQ data 
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Data back to n bit data and 
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Ndata loop 


=1: 


End of for. ipl 


% 


end 


m*pad*Ndata) 


ck (1,m*Ndata* (pad-1) +1 


=Y1p; 


Y2 ((pad-1) *Ndata+1l:pad*Ndata,1) 


end 
end 


hd o£ it 


E 
ih 





% 


N loop 


=1: 


End. .of for b 


% 


.. 


length (ck) ) 


bk=bk (1,1 


Compares rec data string to trans 


% 


¥ 


max (cumsum(err) ) / (m*Ndata*N) 


num2str (BI 


R= 
figure (3) 
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Plot recieved I-Q data with noise 
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axis ( 


grid on 
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